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Summary 



We study the flavor-dependent neutrino spectra formation in the core of 
a supernova (SN) by means of Monte Carlo simulations. Several neutrino 
detectors around the world are able to detect a high-statistics signal from a 
galactic SN. From such a signal one may extract information that severely 
constrains the parameter space for neutrino oscillations. Therefore, reliable 
predictions for flavor-dependent fluxes and spectra are urgently needed. 

In all hydrodynamic simulations the treatment of u^ T and u^ T is rather 
schematic, with the exception of the most recent Garching-Group simulation, 
where the interaction processes were updated partly based on our results. The 
interactions commonly included in traditional simulations are iso-energetic 
scattering on nucleons Nz/ — > z/N, scattering on electrons and positrons 
e ± z/ — > ue ± , and electron-positron pair annihilation e + e~ — > v^v^. In our 
Monte Carlo simulations we vary the static stellar background models in ad- 
dition to systematically switching on and off our set of interaction processes, 
i.e., recoil and weak magnetism in Nz/ — > z/N, scattering on e ± and z/ c /z/ c , 
e + e~ and u c u c annihilation into v^v^ pairs, and neutrino bremsstrahlung 
off nucleons NN — > NNz/z/. As and v^ T sources, NN — > NNz/z/ and 
z/ e z/ e — > v^v^ T- dominate. The latter process has never been studied before 
in the context of SNe and turns out to be always more important than the 
traditional e + e~ annihilation process by a factor of 2-3. For energy transfer, 
the most important reactions are Nz/ — > z/N with recoil, and scattering on 
e ± . Weak magnetism has a very small effect and scattering on z/ c and z/ c is 
negligible. 

The charged current reactions e~p — > z/ e n and e + n — > z/ e p dominate for z/ c 
and z/ e . In comparing our numerical results for all flavors we find the standard 
hierarchy of mean energies (e„ e ) < (e Pc ) < (e^), with, however, very similar 
values for u^ T and z/ e . The luminosities of u^ T and z/ c can differ by up to 
a factor of 2 from L Pc rts L Uc . The Garching Group obtains similar results 
from their self-consistent simulation with the full set of interactions. These 
results are almost orthogonal to the previous standard picture of exactly 
equal luminosities of all flavors and differences in mean energies of up to a 
factor of 2. 

Existing concepts for identifying oscillation effects in a SN neutrino signal 
need to be revised. We present two methods for detecting the earth-matter 
effect that are rather independent of predictions from SN simulations. An 
earth-induced flux difference can be measured by the future IceCube detector 
in Antarctica and a co-detector like Super- or Hyper-Kamiokande. At a single 
detector with high energy resolution the Fourier transform of the inverse- 
energy spectrum can reveal the modulations of the spectrum. Both methods 
are sensitive to the small mixing angle 9i 3 and the neutrino mass hierarchy. 



i 



ii 



Fur meine Eltern 



iii 



iv 



Acknowledgments 



Many people contributed to this dissertation with their advice, support, pa- 
tience, and kindness. I express my gratitude to everyone who helped making 
this work come true. 

Foremost I am indebted to my advisor, Georg Raffelt, not only for propos- 
ing the project and thereby getting me into Astroparticle Physics, but even 
more for his great support and advice in so many ways. I benefited very much 
from our close and very inspiring collaboration. Also for the proofreading and 
uncountable valuable suggestions on the manuscript I am very grateful. 

The numerical part of this work was only possible with the code that 
Hans-Thomas Janka had developed and kindly introduced me to. I enjoyed 
the collaboration with Thomas and appreciate his patience with my ques- 
tions and his exhaustive answers on all kinds of code- and supernova-related 
questions. Thomas and his Garching Supernova Group, i.e., Robert Buras 
and Markus Rampp, contributed very valuable results from their hydrody- 
namic simulations. I also thank Robert for our many long discussions and 
his valuable suggestions on parts of the manuscript. 

With Amol Dighe I had a wonderful collaboration on the neutrino- 
oscillation part. But also for all other questions and discussions I found 
Amol's door wide open. In addition, I thank Amol for improving parts of the 
manuscript by his proofreading and suggestions. Also for discussions with 
Michael KachelrieB, Sergio Pastor, Dimitry Semikoz, and Ricard Tomas I am 
very grateful. 

As important as the support at the institute were all those friends and 
my family who supported and encouraged me throughout all the years, even 
when I was far away. For their love and support in any possible way I owe 
my deepest thanks to my parents as well as my brother and sister. Christoph 
even fitted some proofreading into his busy conference schedule, thank you. 

For your love, support, and encouragement, your patience, your smiles, 
and everything you did for me, words cannot express my gratitude and love 
for you, Meike. 



v 



Parts of this dissertation have already been published as: 



M. T. Keil, G. G. Raffelt, and H. T. Janka, 

"Monte Carlo study of supernova neutrino spectra formation," 

Astrophys. J. in press, 

|astro-ph/0208035| 

R. Buras, H. T. Janka, M. T. Keil, G. G. Raffelt, and M. Rampp, 
"Electron-neutrino pair annihilation: A new source for muon and tau 
neutrinos in supernovae," 



Astrophys. J. 587 (2003) 320 astro-ph/0205006|. 



M. T. Keil, G. G. Raffelt, and H. T. Janka, 
"Supernova Neutrino Spectra Formation," 
Nucl. Phys. B Proc. Suppl. 118 (2003) 506. 

G. G. Raffelt, M. T. Keil, R. Buras, H. T. Janka, and M. Rampp, 
"Supernova neutrinos: Flavor-dependent fluxes and spectra," 
astro-ph/0303226| 



A. S. Dighe, M. T. Keil, and G. G. Raffelt, 

"Detecting the neutrino mass hierarchy with a supernova at IceCube," 
JCAP, in press (2003), j|hep-ph/0303210| 

A. S. Dighe, M. T. Keil, and G. G. Raffelt, 

"Identifying earth matter effects on supernova neutrinos at a single 

detector," 

hep-ph/0304"T50 



vi 



Contents 



1 Introduction! 



2 The Core- Collapse Paradigm! 



2.3 Reincarnation: A Proto-Neutron Star Emerges! 
2_J^_Si;lfiCj2nsis^^ ~ ~~~ 

2.5 Flavor-Dependent Neutrino Emission 



Neutrino Interactions! 

3^^et£rprocessesl . . 



3.2 Neutrino-Nucleon Scattering 

3*3^^^a1^Magnetisml * 

3"^^^remssfaahlung| 

^^_^ajr_^Amrilril^tiQn 

3.6 Scattering on Electrons and Electron Neutrinos! 



IL2 Other Reactions 



Setting the Stage for Neutrino Transport 



^^^^^^^^^^^^^^^^^^^ ^pri^Stars! 



4.2 Proto-Neutron-Star Profiles 



4.3 Optical Depth vs. Thermalization Depth! . 

4.4 Thermalization Depths in Our Stellar Models! 



CJj^ractejj^uij ^j^n^Thermal Neutrino Spectral 



5.1 ( global Parameters! 



5.2 Analytic Fits! 



Relative Importance of Neutral- Current Reactions! 



6.1 Monte Carlo Setup 



6.2 Accretion-Phase Model I 



vii 



6.3 Steep Power Law 



6.4 Shallow Power Law 



IL5 Summary 



Comparison of All Flavorsl 



7.1 Results from Our Monte Carlo Study 



7.2 Previous Literature! 

7.3 Spectral Shapel 

7.4 Summary! . . . . . 



8 Detecting Oscillations of SN Neutrinos! 



Oscillations of SN Neutrinos! 



8.2.1 The Basic Principle] 



Detecting Oscillations With Two Distant Detectors! 



8.2.2 The SN Signal at IccCube 



8.2.3 The Oscillation Signal at Ice Cube 



8.2.4 Super- or Hvper-Kamiokandc and IccCubc 



8.3 Detecting the Earth-Matter Effect at a Single Detector 



£>A Summary 



9 Discussion and Summary 



A Abbreviations! 



R Monte Carlo Code 



B.l General Concept 



B.2 Structure of the Code! 



viii 



Chapter 1 
Introduction 



The physics of neutrino oscillations has received tremendous attention dur- 
ing the past years and today oscillations are thought to be established. A 
large number of experiments has contributed to our present knowledge of 
the mixing schemes and parameters. For several of the mixing parameters 
the allowed ranges are known. Within the next few decades the magnitude 
of the small mixing angle $13 and the question whether the neutrino mass 
hierarchy is normal or inverted will remain open and can be settled only by 
future precision measurements at dedicated long-baseline oscillation exper- 
iments (Barger et al. 2001, Cervera et al. 2000, Freund, Huber, & Lindner 
2001) or the observation of a future galactic supernova (SN). 

Along with the interest in neutrino oscillations developed the branch 
of neutrino astronomy. For the early experiments Raymond Davis Jr. and 
Masatoshi Koshiba received last year's Physics Nobel Prize "for pioneering 
contributions to astrophysics, in particular for the detection of cosmic neu- 
trinos" (Nobel e-Museum, 2002). Current and planned experiments will be 
able to record some 10,000-100,000 neutrinos from a future galactic SN and 
therefore would provide valuable information on SN physics and neutrino 
properties (Barger, Marfatia, & Wood 2001, Minakata et al. 2002). 

There exists a large number of publications addressing the question, what 
the neutrino signal of a galactic SN will tell us about neutrino oscillations. 
With the measurement of such a neutrino signal it would be possible to 
differentiate between existing oscillation scenarios (Chiu & Kuo 2000, Dighe 
& Smirnov 2000, Dutta et al. 2000, Fuller, Haxton, & McLaughlin 1999, 
Lunardini & Smirnov 2001b, 2003, Minakata & Nunokawa 2001, Takahashi 
& Sato 2002). At the same time, an understanding of neutrino oscillations is 
crucial for interpreting a recorded signal. 

In 1987 neutrinos from a SN were detected for the first and only time 
(Bionta et al. 1987, Hirata et al. 1987). This signal from SN1987A in the 
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Large Magellanic Cloud was analyzed taking neutrino oscillations into ac- 
count (Jegerlehner, Neubig, & Raffelt 1996, Kachelriess et al. 2002, Lunardini 
& Smirnov 2001a, Smirnov, Spergel, & Bahcall 1994). However, the neutrino 
detectors that were operational at that time measured only about 20 events 
connected to SN1987A. With such low statistics, determining SN and neu- 
trino parameters was difficult, but some limits could be obtained. 

For the analysis of a SN neutrino signal it would be helpful to have a firm 
prediction for the neutrino spectra emitted by the core of the SN. Oscillation 
effects in a neutrino signal depend on the spectral differences between the 
various flavors. Such spectral differences are obtained by neutrino-transport 
simulations coupled to self consistent hydrodynamic calculations. In order 
to get reliable neutrino spectra the way neutrinos are transported in these 
simulations plays a crucial role. In the past, the transport of was very 
schematic (here and in the following z/ M stands for either muon or tau neu- 
trinos, or anti- neutrinos, unless stated differently). Spectra obtained for 
in these simulations therefore inherit large uncertainties. However, these un- 
certainties barely affect the dynamics of the explosion. From the SN model 
builders point of view, details of spectra formation do therefore not play 
a crucial role. The greatest challenge for numerical models still is to under- 
stand how SNe explode. The most elaborate simulations agree rather well on 
the early stages of the SN mechanism, but do not show explosions. 

According to the commonly accepted "delayed-explosion scenario" a mas- 
sive star with a mass greater than about 8M Q eventually becomes a core- 
collapse SN. At the end of its life, such a star has accumulated an iron core 
of about 1.5 Mq, close to the Chandrasekhar limit. Once the core mass is 
greater than the stability limit it collapses within some hundred milliseconds 
to a proto-neutron star. Due to this collapse a large amount of gravitational 
binding energy is liberated, but particles cannot escape the dense medium. 
A shock wave that formed at the edge of the proto-neutron star travels out- 
wards and plows through the still infalling outer layers of the iron core. The 
disintegration of these infalling iron nuclei consumes the kinetic energy of the 
shock wave and finally causes it to stall after a few 100 ms. Neutrino emission 
drains energy from the proto-neutron star and at the same time heats up the 
region behind the stalled shock wave. This energy transfer causes the shock 
wave to regain velocity and finally blow off the envelope of the star. 

In a crude approximation the proto-neutron star is a blackbody source 
for neutrinos of all flavors. The flavor dependent differences, that are most 
important for oscillations, arise due to the flavor dependence of the neutrino- 
matter interactions. For the case of electron neutrinos (u e ) and electron anti- 
neutrinos (P e ) the dominant reactions are the charged-current interactions 
with nucleons e~p <-> z/ c n and e + n <-> z/ c p. Since there are more neutrons 
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than protons in a proto-neutron star, the interactions are more efficient for 
u c and keep u e in local thermal equilibrium (LTE) up to larger radii than P e . 
Therefore the emerging 9 e are hotter than the u e . The other flavors undergo 
various types of neutral-current interactions that each are of comparable 
importance. The detailed interplay of these reactions in the context of SNe 
has been studied for the first time in the present work. 

Traditionally, the z/^-matter interactions that are implemented in numer- 
ical simulations are iso-energetic scattering with nucleons, scattering on elec- 
trons and positrons, and electron-positron pair annihilation into v^p^ pairs. 
For several years there have been suggestions on improving these rates. Bur- 
rows & Sawyer (1998) and Reddy et al. (1999) calculated the effect of many- 
body correlations that change the neutrino-nucleon interactions. Weak mag- 
netism even causes different rates for and (Horowitz & Li 2000, Horowitz 
2002, Vogel & Beacom 1999). A standard simplification in numerical simula- 
tions was to take the nucleon mass to be infinite for neutrino-nucleon scatter- 
ing. However, it has been known for some time that the influence of nucleon 
recoils in these scattering reactions was not negligible (Janka et al. 1996, 
Raffelt 2001). An important source term for is nucleon bremsstrahlung 
NN — > NNz/^z/^ that even dominates for low energies (Hannestad & Raffelt 
1998, Suzuki 1991,1993, Thompson, Burrows, & Horvath 2000). The tradi- 
tionally included source for in hydrodynamic simulations, e + e~ — > u^u^, 
turns out to be rather unimportant. For the first time we showed that 
electron-neutrino pair annihilation v c v c v^v^ is always a factor of 2-3 
more important than this traditional process. 

In order to study the relative importance of all possible interaction pro- 
cesses and their dependence on a reasonable range of stellar background 
models we have adapted the Monte Carlo code of Janka (1987, 1991) and 
added new microphysics to it. We go beyond the work of Janka & Hillebrandt 
(1989a,b) in that we include the bremsstrahlung process, nucleon recoils and 
weak magnetism, u e u c pair annihilation into v^v^ and scattering of on v e 
and P e . With these extensions we investigate the neutrino transport system- 
atically for a variety of medium profiles that are representative for different 
SN phases. Raffelt (2001) has recently studied the spectra-formation prob- 
lem with the limitation to nucleonic processes (elastic and inelastic scatter- 
ing, recoils, bremsstrahlung), to Maxwell-Boltzmann statistics for the neutri- 
nos, and plane-parallel geometry. Our present study complements this more 
schematic work by including the leptonic processes, Fermi-Dirac statistics, 
and spherical geometry. In addition we apply our Monte Carlo code to the 
transport of u c and u c and thus are able to compare the flavor-dependent 
fluxes and spectra. 

With the complete set of neutrino interactions we find a situation that 
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is almost orthogonal to what has been assumed so far. The standard picture 
concerning SN neutrino fluxes and spectra was an exact equipartition of 
the neutrino luminosity among all flavors and a strong hierarchy of the mean 
energies of the different neutrino flavors [e v ) > (e Pe ) > (e Uc ) with a difference 
of up to a factor of 2 between (e Pc ) and (e„ ). These assumptions in turn lead 
to strong oscillation effects in the predicted signal of a galactic SN. Our 
findings suggest that the difference between (e Pe ) and (e^ ) is about 10% 
or less, but the luminosities can differ by a factor of 2 in either direction, 
depending on the explosion phase. 

This new picture is also supported by all hydrodynamic simulations in- 
cluding Boltzmann solvers for the transport of neutrinos. If these simulations 
use the traditional set of neutrino interactions they find differences between 
(e Pe ) and (e„ M ) around 20-30%. Very recently the Garching group imple- 
mented what we found to be the relevant set of interactions in their two- 
dimensional hydrodynamic simulation (Buras et al. 2003b). Their findings 
are equivalent to the results of our Monte Carlo simulation. 

For studies of neutrino-oscillation effects in SN neutrino spectra these 
findings mean a fundamental change of paradigm, because the main focus of 
these studies are the effects on v c and v^. Available detectors will record only 
u c with high statistics. The basic idea was to identify an oscillation effect 
by comparing characteristics of the observed and the predicted P e spectrum. 
If oscillations were present, one would observe deviations, because the v e 
and Pfj, spectra would mix. In addition to the fact that the predictions have 
large uncertainties, these concepts fail once the new findings are taken into 
account. Therefore different methods are needed. 

According to our findings and the most recent results of self-consistent 
SN simulations, a concept for inferring neutrino-oscillation parameters from 
the signal of a galactic SN should be based on the flux differences between 
v c and rather than differences in mean energies. With a SN rate in our 
galaxy of a few per century the concept needs to involve detectors that will 
be operational for several decades. 

With such methods one might be able to determine oscillation parameters 
that are difficult to obtain by earth-based experiments. If the currently dis- 
cussed "neutrino factories" and "super beams" are built they will be able to 
determine neutrino parameters in great detail (Apollonio et al. 2002, Barger 
et al. 2001, Cervera et al. 2000, Freund, Huber, & Lindner, 2001, Huber et 
al. 2003, Huber, Lindner, & Winter 2002). Today, it is not clear whether such 
experiments will become reality. Both concepts can benefit from one another, 
because detectors being built for earth bases experiments can also serve as 
SN detectors and results obtained from one concept can be used to improve 
the other. 
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In view of our findings, we developed two methods for identifying the earth 
matter effect in a SN neutrino signal that only depend on robust features of 
simulation results rather than detailed predictions. The first method involves 
two neutrino detectors at different sites, such that for one the neutrino path 
goes through the earth and the other sees the signal "from above." The energy 
deposited in one detector differs from that in the other if the earth-matter 
effect is present, due to different original v e and fluxes. Surprisingly, the 
future high-energy neutrino telescope IceCube at the South Pole, that will 
begin taking data from 2005, will yield the most accurate signal (Dighe, 
Keil, & Raffelt 2003a). Due to its unique location, with a co-detector on 
the northern hemisphere, e.g., in Japan, the setup covers most of the sky. 
Alternatively, with a single detector one can identify the modulations of the 
energy spectrum that are induced by the earth-matter effect with the help 
of a Fourier transform (Dighe, Keil, & Raffelt 2003b). This is a difficult task 
with the size and energy resolution of available detectors. However, with a 
detector like Hyper-Kamiokande or a large scintillation detector this is a 
powerful procedure. 

This dissertation is structured as follows. In Chapter |2] we explain the 
current explosion paradigm that we believe is realized in nature. We then 
discuss the neutrino-matter interactions that are relevant inside the proto- 
neutron star in Chapter El We put special emphasis on the u e u e process that 
has never been studied before in the context of SNe. In Chapter 0] we present 
the stellar profiles that we use in our numerical simulations. We introduce 
the concept of a "thermalization depth" and apply it to our stellar profiles 
in order to estimate the relative importance of the interaction processes. 
For describing the spectra that we obtain from our calculations we use the 
characteristic quantities and fits that we explicate in Chapter With all the 
necessary tools in hand we then first describe our numerical results for in 
Chapter |H1 and then move on to the comparison of all flavors in Chapter 
where we also give an overview of the previous literature. As an application 
of our findings we show in Chapter |H] how to identify the earth- matter effect 
in a detected SN neutrino signal. And finally we summarize our findings in 
Chapter H 
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Chapter 2 

The Core-Collapse Paradigm 



Massive stars accumulate iron in their center after several stages of nuclear 
burning. The iron core becomes unstable once it reaches a certain mass limit 
and collapses. According to the so called "delayed-explosion mechanism," 
the collapse is inverted into an explosion and a shock wave travels outwards 
powered by the energy deposition of the neutrino flux. The neutrinos are 
emitted by the proto-neutron star that formed in the center. 

There exist a number of numerical simulations. Even the most elabo- 
rate among these simulations, however, fail to produce explosions. It is an 
unsettled issue what the missing ingredients are. 

With a focus on the neutrino emission, we discuss the various stages of 
the explosion. Flavor dependent differences in the emitted neutrino spectra 
arise due to neutrino-matter interactions. Electron neutrinos and electron 
anti-neutrinos are dominantly produced by charged-current processes. Sev- 
eral neutral current processes govern the transport of the other flavors. 

2.1 The Life of a Star: Balancing Forces 

During its whole life a star has to balance two possibly fatal effects that 
work in opposite directions. On one side, there is the gravitational pull, that 
tries to collapse the star, on the other side the thermal pressure that expands 
the star. The connection of thermal and gravitational energy is given by the 
virial theorem. For the time-averaged kinetic energy (Skin) and gravitational 
energy (Eq) of an atom inside the star the virial theorem is 

(E kin ) = ~(E G ) . (2.1) 

By making the star more compact, (Eg) becomes more negative. Since the 
temperature is proportional to (-Ekin), it rises as the star's radius shrinks. 
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Shrinking the star by heating it up, or vice versa, corresponds to a nega- 
tive heat capacity and stabilizes the star. For a certain stellar radius the 
gravitational pull is balanced by the thermal pressure. 

A star radiates photons and therefore loses energy. By burning just enough 
of its fuel in order to produce the energy radiated away, the star becomes a 
stable object. Initially, this fuel is hydrogen, that produces energy in nuclear- 
fusion reactions to helium. The stronger the gravitational pull, the hotter the 
star's interior and the more energy is radiated away. Therefore more massive 
stars burn fuel at a higher rate. 

Due to the nuclear-fusion reactions more and more helium collects in 
the central region. When hydrogen is depleted at the center, the helium 
eventually reaches high enough density and temperature in order to ignite. 
Due to the strong temperature dependence of nuclear fusion reactions, helium 
burning in the central region is spatially well separated from the hydrogen 
burning that can still continue at the edge of the helium core. The helium 
fusion process produces again heavier elements that can also start fusion 
reactions at a later stage. Since the star has to be denser and hotter to ignite 
the burning of heavier elements it becomes brighter and loses more energy at 
this stage. Additionally, these burning processes liberate less energy in each 
reaction. The reaction rate rises drastically. While hydrogen burning can last 
millions or even billions of years the silicon burning stage usually takes less 
than a day (see e.g. Herant et al. 1997). 

After several stages of nuclear burning the initial composition of helium 
and hydrogen transformed into an onion-shell structure with heavier elements 
up to iron-group nuclei forming the star's core (Figure f2.1j) . For stars with 
a mass greater than about 8 M the core mainly consists of 56 Fe and 56 Ni. 
Due to silicon burning the iron core eventually reaches a mass close to its 
Chandrasekhar limit of 1.2-1.5M . Iron cannot ignite and further fuse to 
heavier elements since it is already the most tightly bound element. 

The Chandrasekhar limit is a stability criterion for compact objects like 
white dwarfs or the iron core of a massive star. Such objects are stabilized 
by electron degeneracy pressure. The following very basic explanation of the 
Chandrasekhar limit can be found in Shapiro & Teukolsky (1983). The grav- 
itational energy of a particular electron at radius R is -Eg —NR' 1 for iV 
gravitating particles. By moving to the center the electron gets lower gravi- 
tational energy and therefore occupies an energetically favored state. As long 
as the electron is non-relativistic its Fermi energy Ep is proportional to the 
Fermi momentum pp squared. Therefore, when the electron moves towards 
the center, Ep oc pp oc N 2 ^ 3 R~ 2 rises faster than the potential energy falls. 
For some radius the total energy Eq + Ep assumes a minimum value and 
the star is in a stable state. When the mass rises and thus the gravitational 
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Fig. 2.1. — Schematic picture of a typical progenitor's onion structure (not 
drawn to scale). 

potential gets deeper, also pp has to rise and more and more electrons be- 
come relativistic. For relativistic electrons, however, Ep oc pp oc A^ 1//3 i? _1 . 
The minimum in the total energy disappears. Depending on its mass being 
above or below the Chandrasekhar limit at that time the star can have two 
distinct destinies. 

• If the mass, M oc N, of the star is below the Chandrasekhar limit 
the total energy is positive. Therefore the star lowers its energy by 
expanding, which in turn lowers pp, electrons become non- relativistic 
and the star is stable again. 

• For a mass above the Chandrasekhar limit the total energy becomes 
more and more negative by lowering the radius and the star starts to 
collapse. 

The number of nucleons determines the mass of the iron core. Up to details 
of the chemical composition of the star, M = Nm nnc \ eon . The Chandrasekhar 
limit is then obtained by solving Ep = Eq for relativistic electrons and yields 
1.5M«. 
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2.2 The Death of a Massive Star: Core 



Collapse 



As soon as the iron core of a SN progenitor reaches its Chandrasekhar limit 
the gravitational pull wins over the internal pressure and the core starts to 
collapse. The rise in density and temperature due to the collapse sets the 
stage for processes speeding up the infall. Electron captures on nuclei and 
protons, 



become more and more efficient and do not only take away degeneracy pres- 
sure, but also drain energy due to the emitted neutrinos that escape the iron 
core unhindered. Another energy sink is photo disintegration of iron nuclei. 
Due to the rising temperature inside the core photons start to disintegrate 
iron nuclei. Within milliseconds the photo disintegration undoes parts of the 
work of the nuclear-fusion processes, and therefore takes away energy and 
thermal pressure. 

As the collapse of the iron core speeds up, two distinct regions develop 
(Figure O left, Figure IP iL. in phase 1): 

• the homologously infalling inner part where the speed is sub-sonic and 
proportional to radius (Bruenn 1985). This region contracts, keeping 
its relative density profile, until it reaches nuclear density where the 
equation of state stiffens drastically. With its enormous inertia the core 
contracts even beyond the new equilibrium and slightly expands again, 
usually referred to as the "core bounce." 

• the super-sonically infalling outer part. Due to the lower sound speed 
at larger radii, the matter assumes super-sonic speeds. The infall cor- 
responds therefore basically to a free fall. 

A shock wave forms at the interface of the two regions due to the slight 
reexpansion of the inner core. The core bounce triggers the explosion. 

The collapse happens on a time scale of some 100 ms (first phase in 
Figure E2J) • Therefore, the outer part of the iron core is still infalling and the 
star's envelope has "not even noticed" that the core collapsed. 



e~(Z,A) 



u e (Z-1,A) 



e p 




2.3. Reincarnation: A Proto-Neutron Star Emerges 
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Fig. 2.2. — Cartoon of a collapsing star. Left: The initial stage, with the 
inner homologous and the outer super-sonically infalling region. Middle: A 
hard core has formed, the shock (white) has bounced and travels outwards. 
The outer layers are still infalling, passing through the shock wave, and then 
dropping at a slower rate. Right: The shock wave has blown off the outer lay- 
ers of the star; the naked proto-neutron star contracts and cools by neutrino 
emission. 

2.3 Reincarnation: A Proto-Neutron Star 
Emerges 

The shock wave moves outwards through the infalling outer layers of the 
iron core and loses energy by disintegrating the nuclei it plows through. The 
nuclear binding energy of 0.1M Q of iron is about 1.7 x 10 51 erg, comparable 
to the explosion energy, i.e., approximately 1% of the released 10 53 erg of 
gravitational binding energy. The dissociation of iron also liberates protons 
that capture electrons, and in turn emit electron neutrinos. In the central 
region of the core, these neutrinos are trapped by the surrounding material. 
Diffusion is slow compared to the dynamic time scale. As shown in Figure l2~3| 
the region of neutrino trapping below R v builds up during collapse. As long as 
the shock wave has not passed R v neutrino trapping is very efficient, because 
there are still iron nuclei around and due to coherence effects the cross section 
for neutrino scattering on nuclei is large compared to that for scattering on 
nucleons. 

Then we enter the second phase of Figure 12.31 The shock passes R v and 
disintegrates the still infalling iron nuclei in the vicinity of R v . Suddenly, there 
are no longer iron nuclei around. Therefore, instead of neutrinos scattering 
on nuclei, only scattering on free nucleons takes place with a much lower 
cross section. At the same time, the free protons capture electrons and emit 
electron neutrinos. These neutrinos in the vicinity of the shock can escape 
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Fig. 2.3.— Schematic plot of the time evolution for various radial zones 
within the radius of the original iron core. The evolution is divided in four 
phases: 1. Collapse of the iron core, 2. prompt z/ e burst, 3. accretion phase, 
and 4. Kelvin-Helmholtz cooling of the proto-neutron star. As long dashes the 
radius of the iron core (RFe) is shown, short dashes represent the interface to 
the inner core (-Ri C ), and dots to the radius below that neutrinos are trapped 
(R v ). The solid line, starting shortly after 0.1 s shows the position of the 
shock (-R s hock)- Adapted from Janka (1993). 

freely and cause a sudden peak in the electron-neutrino flux. The luminosity 
briefly rises to around 10 54 erg s _1 , referred to as the prompt neutrino burst. 
This also weakens the shock. 

Lepton number is approximately conserved inside the core because neu- 
trinos are trapped in the dense matter mainly due to neutrino-nucleon scat- 
tering. The processes of Equations ()2.2j) cause a (3 equilibrium to develop and 
a Fermi sea of v c builds up. Most of the gravitational binding energy is now 
stored in the Fermi seas of electrons and electron neutrinos. 

During the third stage in Figure 12.31 at a radius of a few 100 km, only 
about a few 100 ms after bounce, the shock wave stalls due to its energy 
loss caused by the disintegration of iron nuclei and neutrino emission. Still 
infalling matter from outer layers of the iron core passes through the standing 
shock onto the newly formed proto-neutron star, also shown in the middle 
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panel of Figure l2~2l The edge of the proto-neutron star is at i? ic in Figure l2~3l 
what used to be the interface between the homologously and super-sonically 
infalling regions of the iron core. 

Contraction of the proto-neutron star and the accretion of matter still 
cause the temperature to rise and the neutrino luminosity increases. In the 
region between R u and the standing shock wave convection takes place. The 
neutrinos cool the surface of the proto-neutron star. In the last stage of 
Figure 12.31 the stalled shock wave is revived by the energy that a fraction 
of the outgoing neutrinos deposit in the region just behind the shock. If 
only a few per cent of the liberated 10 53 erg gravitational binding energy are 
absorbed by the matter behind the standing shock wave, the shock will move 
outwards again. Finally, this powerful explosion blows off the whole envelope 
of the star. The naked proto-neutron star remains in the center and cools by 
neutrino emission. The Kelvin-Helmholtz-cooling phase is displayed as the 
last phase in Figure 12.31 and the right panel of Figure 12.21 

Due to the initially stalled and afterwards revived shock wave, this mech- 
anism is referred to as a delayed explosion. It was first suggested by Bethe 
and Wilson (1985). Over the years more and more refined numerical simu- 
lations have been developed. The most elaborate among them agree in their 
result that no explosions are obtained (Buras et al. 2003a, Liebendorfer et 
al. 2001, Mezzacappa et al. 2001, Rampp & Janka 2000, Thompson, Burrows, 
& Pinto 2002). It has been stressed by Buras et al. (2003a) that some crucial 
piece in the game might be missing. 

Even if the standard picture is correct, there are also cases when no 
neutron star is born by the SN explosion. If the progenitor star has a large 
mass, say greater than 20 M , the evolution of the proto-neutron star can 
find a sudden end. Due to accretion the mass of the proto-neutron star can 
exceed its Chandrasekhar limit. The proto-neutron star then collapses and 
ends up as a black hole, terminating neutrino emission almost at an instant. 
The details depend on the nuclear equation of state and matter accretion, 
issues that are not yet settled. A discussion can be found in Beacom, Boyd, 
& Mezzacappa (2001) and references therein. 

2.4 Self-Consistent Simulations 

There have been numerical simulations of SNe showing delayed explosions 
for almost two decades now (Bethe & Wilson 1985, Wilson 1985). These pi- 
oneering works obtained robust explosions. With increasing computer power 
and refined input physics the previous explosions could not be confirmed 
by other groups. Up to now, only the Livermore group obtains robust explo- 
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sions (Totani et al. 1998). However, they have to include convection inside the 
proto-neutron star just below the neutrino spheres by a parameterization of 
neutron-finger instabilities in order to get explosions (Wilson & Mayle 1988, 
1993). This treatment of convection has the effect of enhancing the early 
neutrino luminosities, which in turn causes more energy to be deposited be- 
hind the shock. By neutron-finger convection the neutron rich outer part of 
the proto-neutron star mixes with the inner part in analogy to salt-finger 
convection in hot salt water and cold fresh water. However it is controversial 
whether convection is realized in an actual SN (Burrows 1987, Bruenn, Mez- 
zacappa, & Dineva 1995, Keil, Janka, & Miiller 1996, Mezzacappa et al. 1998, 
Pons et al. 1999). In their two-dimensional simulations, Buras et al. (2003b) 
find convection of Ledoux type. 

All groups performing computer simulations of the hydrodynamics of SNe 
with a state-of-the-art treatment of the neutrino transport, i.e., Boltzmann 
solvers, do not get explosions at all (Liebendorfer et al. 2001, Mezzacappa et 
al. 2001, Rampp & Janka 2000, Thompson, Burrows, & Pinto 2002). Even 
the most elaborate SN simulations with multi-dimensional hydrodynamics 
and all relevant microphysics in the neutrino interactions included fail to 
explode (Buras et al. 2003b). The reason for the models not to explode is 
obscure. 

The delayed-explosion mechanism described in the previous section is in 
principle found by numerical simulations. The collapse sets in, a shock wave 
forms, bounces and moves outwards. On its way it loses energy mainly due 
to the disintegration of nuclei and finally stalls. The simulations also show 
the appearance of a cooling region above the neutrino sphere where energy 
is lost due to neutrino radiation. Behind the shock these neutrinos heat the 
matter, as described above. In the simulations, however, the heating is not 
sufficient for reviving the shock. Instead it falls back and no explosion takes 
place. Quite naturally, the main focus of model builders is to find the missing 
input physics in order to obtain an explosion. 

In spite of the numerical failure of obtaining explosions, there are good 
reasons to believe in the delayed-explosion scenario. Observational hints sug- 
gest that the core collapse paradigm should, at least to some extent, be 
realized in nature. We know that SNe occur after stars have reached their 
final stage of nuclear burning. For example, SN1987A was a blue supergiant 
before it exploded. There were neutrinos observed from SN1987A (Bionta et 
al. 1987, Hirata et al. 1987) about three hours before the first photons were 
seen, which is a clear evidence for a core collapse. Although there has not 
been a neutron star identified at the site of SN1987A, for other SN remnants 
neutron stars have been observed. While the details of the explosion mech- 
anism are still under debate, basic features appear to be robust. If a core 
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Fig. 2.4. — Numerical results of the first 2 s after core bounce taken from 
Totani et al. (1998). Top left: Only the mean energy of u e is affected by the 
prompt burst. Top right: The evolution of mean energies shows the traditional 
hierarchy. Bottom left: The prompt u e burst lasting for a few milliseconds. 
Bottom right: Time evolution of the neutrino luminosities, scaled on both 
axes as compared to the left panel. 



collapse occurs the rise of neutrino luminosity and also the prompt neutrino 
burst are unavoidable as long as the collapse does not continue to form a 
black hole immediately An ongoing neutrino emission for several seconds by 
the proto-neutron star is also not doubted. Throughout this thesis we will 
assume that the delayed-explosion mechanism is realized in nature. 

All numerical simulations agree on this qualitative behavior. An example 
that we took from Totani et al. (1998) is given in Figure E3I Neutrino mean 
energies and luminosities are shown as functions of time post bounce. Solid 
lines represent u e , dashed u e , and dotted v^. In the bottom-left panel we show 
a prominent feature, namely the prompt v e burst right after the bounce. The 
mean energy of u e in the top-left panel shows a short rise during the prompt 
burst. The left panels show the mean energies and luminosities during the 
first 60 ms and the right panels during the first 2 s after core bounce. Right 
after the prompt u e burst the accretion phase takes over, visible in the rise and 
fall of luminosities. After the end of the accretion phase deleptonization and 
cooling of the proto-neutron star powers the neutrino emission. The cooling 
phase starts at 0.5 s. How the flavor-dependent differences arise is the topic 
of the next section. 
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2.5 Flavor-Dependent Neutrino Emission 

As soon as the collapse starts the neutrino luminosity increases. In this early 
stage only electron neutrinos are produced in large amounts due to electron 
captures on free protons. The protons originate from photo-disintegrated iron 
nuclei. Since these processes happen most frequently in the dense inner core 
the major part of the released electron neutrinos is trapped. This trapping is 
the reason for the Fermi sea of u e to develop. When the shock wave passes the 
region where the star becomes transparent to neutrinos, i.e., the "neutrino 
sphere," the trapped v e can escape. Additionally, the infalling iron nuclei are 
disintegrated and the electron capture rate in the vicinity of the neutrino 
sphere goes up. Altogether, this results in the prompt v e burst (left panels 
Figure EPj) . 

Neutrinos below the neutrino sphere are still trapped. During the accre- 
tion phase, lasting for a few 100 ms, infalling matter deposits energy on the 
surface of the proto-neutron star, increasing the neutrino luminosities (bot- 
tom right panel Figure|23J). Afterwards, the proto-neutron star still gains en- 
ergy by contraction. The Fermi seas of electrons and z/ e have a huge amount 
of energy stored. The star releases this energy over a timescale of several 
seconds. For the calculation of the Livermore group (Totani et al. 1998) we 
plotted the deleptonization in Figure 1231 The lower panel shows the number 
of leptons per nucleon (Yl) as a function of radius. From top to bottom the 
lines correspond to 0.5, 1, 3, 5, 10, and 15 s after bounce. Accordingly, the 
top panel shows the temperature profile at these times. The maximum moves 
to lower radii as time increases. At the same time the innermost part heats 
up. 

Even though neutrinos are trapped below their neutrino spheres, due 
to the steep gradient in Y"l they diffuse outwards. Since electrons are very 
degenerate in this region, even more degenerate than v ei neutrinos will scat- 
ter towards regions of lower electron degeneracy. Additionally, when scat- 
tering on electrons, neutrinos will always lose energy, since the electron has 
to end up with an energy above the surface of the electrons' Fermi sea. By 
this down scattering the electron-neutrino- and electron-degeneracy energy 
is mainly released as heat. Therefore the temperature maximum follows the 
deleptonization (see e.g. Raffelt 1996). 

The reason why neutrinos are trapped is their frequent interaction with 
the surrounding matter via various interaction processes. The medium basi- 
cally consists of neutrons, protons, electrons, positrons, and of course neutri- 
nos. Due to the presence of electrons and positrons, z/ e and v e can undergo 
charged current reactions that are absent for all other flavors. Qualitatively 
the transport of v e and v e is therefore different from the transport of and 
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Fig. 2.5. — Deleptonization of a proto-neutron star as obtained by the Liver- 
more group (Totani et al. 1998). Top: We show the temperature as a function 
of radius at 0.5, 1, 3, 5, 10, and 15 s after bounce. With increasing time the 
peak moves to lower radii. Bottom: For the same times we give the lepton 
number per nucleon. As time increases lepton number gets lower. 

v T . To be exact, there is also a small population of muons and anti-muons 
present. In central regions, where the temperature is maximal, the muon 
mass is twice the temperature. Since for a non-degenerate particle the mean 
energy is roughly three times the temperature , (e) ~ 3 T, there is a thermal 
distribution of muons present in the inner core. However, in the vicinity of 
the neutrino sphere of and z/ M , where muons would affect the emerging 
neutrino spectra, the temperature is about 1/10 of the muon mass. Thus 
muons are strongly suppressed and do not affect the neutrino emission. 

Schematically, the emission of z/ e and z/ e is displayed in Figure 12.61 In the 
central region of the proto-neutron star, corresponding to the black shaded 
area on the left hand side of Figure 12.61 /5-processes keep z/ e and z/ e in LTE. 
Due to the decreasing density this reaction becomes inefficient at some ra- 
dius, defining the neutrino sphere. From that radius neutrinos start streaming 
freely. Since the /3-cross section depends on neutrino energy also the radial 
position of the neutrino sphere is energy dependent. Therefore it is not pos- 
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Free streaming 




Fig. 2.6. — Schematic picture of the u e and u e spheres in a proto-neutron 



Fig. 2.7. — Same as Figure EE] but for v^. 

sible to define a unique neutrino-sphere radius. Depending on their energy 
neutrinos decouple from the proto-neutron star at different radii and there- 
fore different local temperatures. For this reason the spectrum of neutrinos 
leaving the star does not represent a thermal distribution with a tempera- 
ture corresponding to the medium temperature at the neutrino sphere that 
is, anyway, not well defined. 

The energy dependence of the interaction processes of v e and v e is exactly 
the same. The proto-neutron star, however, contains more neutrons than 
protons, which leads to a higher absorption rate for u e than for v e . Therefore, 
the u e sphere lies at a larger radius than the u e sphere. Since this is true for 
all neutrino energies the mean energy of emitted v e is always less than the 
mean energy of u e . 

For rough estimates, the neutrino sphere can be considered a blackbody 
radiating neutrinos. The surface area is thus proportional to the number 
of emitted neutrinos. Again, the energy dependence of the neutrino-sphere 
radius makes a more reliable estimate very difficult. 

In Figure 12*771 we show a schematic view of the transport. Comparing 
Figures 12.71 and 12.61 tells us that the transport of is more complicated. 
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There are no charged current interactions of with the medium. Instead, 
there is a variety of neutral current interactions that qualitatively differ from 
one another and, at the same time, are of comparable importance for the 
transport process. Naively, one would expect that decouple deeper inside 
the star than P e , because the neutral-current interaction rates are lower. 
Then, the mean energy of emerging would be higher compared to u c and 
the flux lower. In reality, the situation is more subtle. 

In the inner-most region the are kept in LTE by energy exchanging 
scattering processes and the following pair processes: 

• Bremsstrahlung NN <-> NNz/ M z/ M 

• Neutrino-pair annihilation v c v c <-> v^v^ 

• Electron-positron-pair annihilation e + e~ <-> v^v^ 

The radius where these reactions become inefficient defines the number 
sphere. Outside this sphere no particle creation or annihilation takes place 
and can no longer be in LTE. 

Due to the scattering reactions 

• Nz/ M -> Nz/^ 

• e ± u fl -> e ± z/ M 

there is still energy exchange with the medium. However, the two processes 
are qualitatively very different. Scattering on e^ 1 is less frequent since there are 
less e ± than nucleons. On the other hand the amount of energy exchanged in 
each interaction with e 1 * 1 is very large compared to the small recoil of nucleons. 
At the radius, where scattering on e ± freezes out lies the energy sphere. A 
diffusive regime starts, where neutrinos only scatter on nucleons and therefore 
exchange little energy in each reaction. This regime is terminated by the 
transport sphere, defined by the radius at which also scattering on nucleons 
becomes ineffective and the start streaming freely. 

Even if we assume that nucleon scattering happens iso-energetically, i.e., 
no energy can be exchanged outside the thermalization sphere, the mean 
energy of neutrinos escaping the star is only about 50-60% of those close 
to the thermalization sphere. Due to its dependence on the square of neu- 
trino energy the nucleon scattering cross section has a filter effect, because it 
tends to scatter high energy neutrinos more frequently (Raffelt 2001). If one 
were to specify a neutrino sphere according to the crude approximation of 
a blackbody radiating neutrinos, it would coincide with the number sphere. 
Neutrinos leave LTE at the number sphere. The position of the number sphere 



20 



Chapter 2. The Core-Collapse Paradigm 



determines the flux, because neutrino creation is not effective beyond that 
radius. The flux that passes the number sphere is conserved. On the other 
hand, the mean energy of in this area is still significantly lowered due 
to scattering processes before the leave the star. The mean energy of 
emerging the star is usually found to be larger than that of z/ e . Simulations 
that include the interactions only very approximate obtain a large hierar- 
chy in the mean energies (Figure I23J. The systematic study of all interaction 
processes is a major part of this work. Our findings show that the transport 
of is more subtle than previously assumed. Once all interaction processes 
are taken into account, the mean energies of and z/ e are very similar and 
might even cross over. The number fluxes can differ by large amounts. 



Chapter 3 

Neutrino Interactions 



We discuss neutrino interactions with the stellar medium that we assume 
to consist of protons, neutrons, electrons, positrons, and neutrinos of all fla- 
vors. The most important reactions for electron neutrinos and electron anti- 
neutrinos are charged current processes on nucleons (/9-processes). All other 
flavors undergo neutral-current interactions only. 

The neutral-current processes fall into three qualitatively different groups. 
True LTE is only obtained by the pair creation and annihilation processes, 
i.e., bremsstrahlung, the u e u e , and e + e~ pair processes. Scattering on e 1 * 1 
and z/ e , v e exchanges large amounts of energy in a single reaction. Neutrino- 
nucleon scattering qualitatively differs from the other scattering reactions in 
being more frequent but at the same time exchanging only small amounts of 
energy in each reaction. Going beyond the traditional iso-energetic nucleon 
scattering we introduce recoil and weak magnetism. 

For the first time, we have investigated the z/ e z/ e pair process in detail and 
were able to show that it is far more important than the traditional e + e~ pro- 
cess. By crossing symmetry the v e v e pair process is related to v^Vc scattering. 
This scattering process, however, is negligible compared to e 1 * 1 scattering. We 
have presented these results in a similar form in our publication: 
R. Buras, H. T. Janka, M. T. Keil, G. G. Raffelt, and M. Rampp, "Electron- 
neutrino pair annihilation: A new source for muon and tau neutrinos in su- 
pernovae," Astrophys. J. 587 (2003) 320. 

3.1 Beta-Processes 

For electron neutrinos, the most important neutrino-matter interactions are 
the /^-processes z/ e n <-> e~p and z/ c p e + p. In principle, also can un- 
dergo charged current reactions. The central region of the SN core certainly 
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contains muons, but their presence is commonly neglected in hydrodynamic 
simulations. In the vicinity of the neutrino sphere the temperature is too 
low for a thermal distribution of muons. For our simulations it is therefore 
justified to neglect the charged current reactions of v^. 

The squared spin-summed matrix element for neutrino absorption by neu- 
trons is given by (Yueh & Buchler 1976b) 

J2\M\ 2 = 32G 2 \(a + l) 2 (k 2 -h)(h-k 3 ) 



spins 



+ (a - l) 2 (k 2 ■ k 3 )(k A ■ ki) + (a 2 - 1) m n m p (ki ■ k 3 ) 



(3.1) 



The momenta k\, k 2 , k 3 , and k± are assigned to the corresponding particles as 
shown in Figure I3~T| Gp is the Fermi coupling constant, a = Ca/Cv = 1.26, 
and m n and m p the masses of the neutron and proton, respectively. Together 
with phase-space blocking and the phase-space integrations, Equation (jH.lj) 
yields the absorption rate or equivalently the inverse of the mean free path 
(mfp) 

x n n f d 3 k 2 d 3 k 4 d 3 k 3 2 



(2vr) 5 2ei J 2e 2 2e 4 2e 



3 



spins 



x/ n (e 2 ) (1 - / p (e 4 )) (1 - / e (e 3 )) 5\h + k 2 - k 3 - fc 4 ) , (3.2) 

where n n is the neutron density, ej the energy of particle j, and f n {^i)i / p (q), 
and /e(es) are the occupation numbers of neutrons, protons, and electrons, 
respectively. 

For anti-neutrinos we have to interchange neutrons and protons and re- 
place electrons by positrons. For the reverse rate one has to replace blocking 
factors and the occupation numbers appropriately. 




Fig. 3.1. — Feynman diagram for u e absorption. 
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Since there are more neutrons than protons present in the SN core, the (5- 
rates for u e are always higher than for v e . All interaction processes presented 
in the remaining sections of this chapter are neutral-current interactions. 
Due to the smaller coupling constants they are sub-dominant for electron 
neutrinos. The following sections are more relevant to z/ M , because charged 
current reactions do not contribute to the spectra formation of v^. (If not 
stated differently, always refers to v^ T and v^ T .) 

3.2 Neutrino-Nucleon Scattering 

The most frequent neutral-current process is neutrino-nucleon scattering and 
therefore the dominant opacity source for v^. For a neutrino with initial 
energy e\ and final energy €2, the differential cross section is given by 



with 00 = 6\— £2, k the modulus of the momentum transfer to the medium, and 
the scattering angle. S(u, k) represents the dynamical structure function 
that parameterizes the response of the nuclear medium. 

The dynamical structure function includes all effects that change the be- 
havior of the nucleon in the medium, like fluctuations, correlations, and de- 
generacy effects. The full dynamic structure function covers the whole (u, k) 
plane. Therefore it does not only describe neutrino scattering on nucleons 
but also neutrino bremsstrahlung off nucleons. Both processes are related 
by crossing symmetry as shown in Figure 13 .'2\ where the grey shaded blob 
represents the nuclear medium. However, the full structure function is un- 
known. Common approximations treat the space-like (uj 2 < k 2 ) and time-like 
(oj 2 > k 2 ) domains on a different footing. 



Fig. 3.2. — Left: Diagram of neutrino nucleon scattering. The grey shaded 
blob represents the nuclear medium. Right: Bremsstrahlung is obtained by 
crossing one neutrino leg. 




2 



(3.3) 




24 



Chapter 3. Neutrino Interactions 



The common simplification of the neutrino-nucleon scattering cross sec- 
tion in traditional SN simulations was to assume that the reaction is iso- 
energetic, i.e., the nucleons have an "infinite mass" and therefore cannot 
absorb energy but arbitrary amounts of momentum. In this approximation 
the structure function is 

Sno-recoil^, k) = 2tt5(u) . (3.4) 

The only multi-particle effect that is usually taken into account is final-state 
blocking due to the nucleon degeneracy. 

More elaborate recent approximations include correlated nucleons and 
a finite nucleon mass. (Burrows and Sawyer 1998, Reddy, Prakash, and 
Lattimer 1998, Reddy et al. 1999). These results were calculated in the 
random-phase approximation, that includes nucleon-spin correlations but not 
nucleon-spin fluctuations. Therefore, bremsstrahlung is not allowed, because 
neutrinos couple to the nucleon spin and the emission of a neutrino pair re- 
quires a spin flip. In the recent simulations by the Garching group (Rampp 
and Janka 2002), the correlation effects were taken into account together 
with recoil. 

Some different aspects of multi-particle effects are taken into account 
by nucleon excitations. The diagrammatic structure is given in Figure 13.31 
Raffelt (2001) developed a schematic structure function in this spirit for 
bremsstrahlung. By crossing a neutrino leg one obtains inelastic scattering, 
that allows for an energy transfer unrelated to nucleon recoil. Even though, 
this energy transfer is of the same size as nucleon recoil, the quadratic den- 
sity dependence of the inelastic scattering process yields lower rates in more 
dilute regions, where recoil is still effective. In a detailed study Raffelt (2001) 
shows that once recoil is included in a SN simulation the inelastic scattering 
contribution can be neglected. For the implementation of z/ M -nucleon scatter- 
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Fig. 3.3. — Left: Diagram of neutrino bremsstrahlung off nucleons. This pro- 
cess is subdominant for energy exchange compared to recoil (Raffelt 2001). 
Right: The "crossed process" corresponds to inelastic neutrino-nucleon scat- 
tering. 
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ing Raffelt (2001) provides a structure function for recoiling nucleons that 
we adapted to our code. 

We do not distinguish between protons and neutrons. Since for non- 
relativistic nucleons the scattering cross section is proportional to C v + 3C A , 
the vector current (CV = — \ for neutrons and \ — 2sin 2 6 l w for protons) 
is small compared to the axial component, where we use \Ca\ = 1.26/2. 
Neglecting the vector part simplifies the calculations significantly, since oth- 
erwise, there are different structure functions for the axial current, the vector 
current, and the mixed term. 

Ignoring nucleon degeneracy effects, the structure function that incorpo- 
rates nucleon recoils is (Raffelt 2001) 

W^, k) = ^exp [~^] , (3-5) 

with Uk = k 2 /2m, m the nucleon mass, and T the medium temperature. 

Multiplying Equation ()3.3j) with the density of nucleons and ignoring 
phase-space blocking of the essentially non-degenerate nucleons yields the 
differential rates that can be integrated for obtaining the required energy 
and angular differential rates. In the case of recoil the numerical integra- 
tions are rather tricky because Equation (j3.3J) is strongly forward peaked, cf. 
Figure E3J In our code we employ the "rejection method" for obtaining the 
integrated rates (Press et al. 1992). 




Fig. 3.4. — The differential interaction rate for i/N scattering according to 
Equations ()3.3|) and ()3.5|) for initial energy e\ = 10 MeV as a function of the 
outgoing energy e 2 and scattering angle /x = cos 9. 
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In the case of v e and u e , scattering on nucleons is less important than the 
/3-processes. Therefore we left the original implementation of our code in the 
electron-neutrino sector unchanged. The implementation follows Tubbs & 
Schramm (1975) and can be found in Janka (1987). It basically corresponds 
to using the structure function given in Equation (|3.4|) . 



3.3 Weak Magnetism 

Weak magnetism arises as a correction to nucleon scattering due to our 
schematic treatment of the dynamical structure function. The complete struc- 
ture function would, of course, contain weak magnetism. The effect is absent 
if we use the "traditional" approximation of infinite nucleon mass. 

As mentioned earlier, there are essentially no muons or taus present in the 
vicinity of the neutrino sphere. Therefore, and are usually treated in the 
same way. However, the neutrino-nucleon cross section is different for neutri- 
nos and anti-neutrinos once weak magnetism is included. The contributions 
to the cross section from the anomalous magnetic moment of nucleons and 
the interference term of the axial and vector current have different signs for 
neutrinos and anti-neutrinos due to parity violation in the standard model 
(Horowitz 2002). The interaction rate is higher for neutrinos than for anti- 
neutrinos. 

Neutrinos in the relevant regime of the proto-neutron star have energies 
on the order of several 10 MeV. To first order in e/m, where e is the neutrino 
energy and m the nucleon mass, the correction factor to the neutrino-nucleon 
cross section due to weak magnetism is (Horowitz 2002) 

AC A (C V + F 2 ) k 

C 2 A (3-cos6) m' 1 J 

where k is the momentum transfer, F 2 a structure function parameterizing 
the structure of nucleons. For neutrino-proton scattering F 2 = |(/i p — // n ) — 
2 sin 2 (6 l w )/ip « 1.019, where // PjI1 are the magnetic moments of neutron and 
proton, respectively, and 8 W is the weak mixing angle. For neutrino-neutron 
scattering, the indices n and p have to be interchanged yielding F 2 = —0.963. 
The upper sign in Equation ()3.6|) is to be used for neutrinos and the lower 
sign for anti-neutrinos. 

In order to be consistent with Equation ()3.3|) we keep terms cx C\ and 
neglect terms oc C v . In addition, we substituted eicos^, taken at the rest 
frame of the nucleon, by our momentum transfer k. This is correct for for- 
ward and backward scattering but only an approximation for other angles. 
It does not change the angular dependence by much, but in order to keep 



3.4. Bremsstrahlung 



27 



our structure-function prescription of the cross section, we cannot multiply 
factors that depend on k. 

The first-order correction factor has the disadvantage that it becomes 
negative for large momentum transfer. In order to avoid numerical problems 
we use 

( 1± iC A (Cy + F 2 ) kV 

\ C 2 A (3-cos6) 2m) { ' 

instead, which corresponds to Equation (|3.6|) up to first order in k/m. For 
neutrino momenta in the region of the neutrino sphere this approximation 
induces an error of less than 1%. 



3.4 Bremsstrahlung 

Neutrino bremsstrahlung off nucleons inside the stellar medium is related 
to neutrino nucleon scattering by crossing symmetry (Section 13.2)1 . Raffelt 
(2001) showed that the emerging neutrino spectra from a SN do not sensi- 
tively depend on the detailed rate, a fact that we confirm in our simulations. 
A schematic approach is therefore well justified. Since bremsstrahlung only 
plays a crucial role for we do not consider u e and v e . We adopt the ap- 
proach given in Raffelt (2001). The rate for the absorption of a by inverse 
bremsstrahlung NNz/ M z/ /t — > NN is given by 

A" 1 = ^ n B ± / m 24 e eS(cj) , (3.8) 



2 2e J 2e(2vr) 3 

where uj = e + e, ti-q is the nucleon density, and \Ca\ = 1.26/2 as in the 
scattering case. The over-barred quantities belong to the that is absorbed 
together with the primary v^. 

The energy-differential rate for emission of bremsstrahlung is obtained by 
adjusting the phase-space blocking: 

Hr) r 2 Ci 2 p 2 r rl 3 k 

t - -f" - h [1 " /(e)1 / Wik 11 - /(£_)1 24 "" SM (M) 

The heuristic ansatz for S(u) given by Raffelt (2001) has the form of a 
Lorentzian 

g ("> = i + 7 m - r = f)t" Br 2 - < 3 ' 10 ' 

lj 2 + T 2 l+exp(-w/T) m 2 JmT + m 2 



with the constant a n ~ (2m/m 7r ) 2 /47r, where m is the nucleon mass and 
the pion mass. For carrying out the integration in Equations (|3.8)1 and ()3.9|) 
the r 2 in the denominator of S can be neglected. A very detailed analysis of 
this bremsstrahlung rate can be found in Raffelt (2001). 
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3.5 Pair Annihilation 

We now turn to the leptonic source reactions for v^v^ pairs, e + e~ — > v^v^ 
and z/ e z/ e — > v^v^. Traditionally, only e + e~ pair annihilation was implemented 
in SN simulations. As one part of this dissertation we have for the first time 
investigated the importance of z/ e z/ e pair annihilation that turns out to be a 
factor of 2-3 more important (Buras, Janka, Keil, Raffelt, & Rampp 2003a). 

The matrix elements for both processes are identical up to coupling con- 
stants while the phase-space integrations only differ by the chemical poten- 
tials. After summing over all spins and neglecting the rest masses, the squared 
matrix element is 

J2 \M\ 2 = 8 Gl[{C v + C A fu 2 + (C v - C A )H 2 ] (3.11) 

spins 

with the Mandelstam variables t = —2k\ ■ k% and u = —2k\ ■ k^. The momenta 
are assigned to the particles as indicated in Figure 13*31 The weak interaction 
constants for e + e~ annihilation are 

CV = -i + 2sin 2 fl w , C a = -\ (3.12) 
while for z/ e z/ e annihilation they are 

C A = C V = ±. (3.13) 

For the interaction rates we have to perform the phase-space integrations, 
using blocking factors for the final states and occupation numbers for initial- 
state particles (Hannestad & Madsen 1995, Yueh & Buchler 1976a). Three 
integrations remain that can not be carried out analytically. 

In order to perform the phase-space integrals we have to specify chem- 
ical potentials for all interacting particles. Mu- and tau-leptons are almost 




Fig. 3.5. — Pair annihilation processes producing v^v^ pairs. 
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absent in proto-neutron star atmospheres, therefore, chemical potentials of 
the corresponding neutrinos can only arise due to differences in the neutrino 
and anti-neutrino interactions, i.e., weak magnetism. This chemical potential 
builds up dynamically, because anti-neutrinos escape more easily than neu- 
trinos. For the e + e~ reactions the local value of /x e can be obtained from p, 
T, and Y e by inverting 

n e -(/i e ) -ne+Hjj =Y (3 14) 

^■baryons 

where n e -(/i e ) and n e +(— /i e ) are the number densities of electrons and 
positrons, respectively. The baryon number density is 

^-baryons 7Z 7TT ' 7T" ■ (3.15) 

m n (l - Y c ) + m p Y e 
For u e and v e the chemical potential is obtained by the relation 

Pu c = Pe + Hp- Pn, (3.16) 

with the chemical potentials /i p and /x n of protons and neutrons, respectively. 

The reason why the z/ e z/ e pair process is an important source for v^v^ 
pairs is that below the number sphere, v e and V e are basically a part of 
the stellar medium like e 1 * 1 . The energy sphere of v ^ lies always deeper inside 
the star than the u e and u e spheres (Section I4.4j) . Thus u e and u e are in LTE 
and are part of the medium as far as the transport of is concerned. 

Until very recently e + e~ annihilation was considered to be the most im- 
portant leptonic source for v^. In most numerical simulations it is still the 
only included creation mechanism for the flavors. Comparing the rates for 
SN conditions it is obvious that v e v e pairs are more important. In Figure I3~6l 
we show both leptonic pair rates as a function of their respective degeneracy 
parameters. The chemical potential of was set to zero and u e , i> e , and e 1 * 1 
are assumed to be in LTE. In the region below the neutrino spheres the de- 
generacy parameters obey i] Uc = r] e + r] p — r] n < n e and therefore the neutrino 
pair process is always more important than e + e~ annihilation. 

Another qualitative difference of both rates becomes apparent by looking 
at the energy-differential production rate. In the case of vanishing degeneracy 
of e 1 * 1 and z/ e , v e the differential production rate d 2 n/dedt is equal for both 
and as shown in Figure I3~71 In general, however, electrons and electron 
neutrinos will be significantly degenerate. For the extreme case rj e = 7] Va = 10 
we show the differential production rates in Figure 13.81 In the upper panel 
both and creation rates are very similar for e + e~ annihilation. This 
represents the fact that the prefactors u 2 and t 2 in Equation 1)3.11)1 are similar, 
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Fig. 3.6. — Pair production rates by the process v e v e — > v^v^ as a function 
of r/ Uc (upper line) and e + e~ — > v^v^ as a function of rj e (lower line). We used 
T = 12 MeV and rj Vl = 0. 



namely (CV + Ca) 2 ~ 0.54 2 and [C\ — Ca) 2 ~ 0.46 2 . Comparing the rates for 
and u^, which corresponds to exchanging u and t, then has no big effect. 
For the u e u e process the situation is different. In this case (Cy + Ca) 2 = 1 
and (CV — Ca) 2 = and only u 2 contributes. Interchanging z/ M with 
corresponds to replacing u 2 by t 2 . Therefore, the kinematics for particles and 
anti-particles are different. Since in the case of n Uc ^ also the distribution 
functions of particles and anti-particles are different this leads to a difference 
in the energy dependence of the and production rate. The lower panel 
of Figure 13.81 displays both rates. The numbers of produced and are 
equal, of course. 

The reason for the average energy of to be larger than of can be 
understood by looking at the reaction v e v e — > v^v^ in the center of momentum 
(CM) frame. The differential cross section is 

d(T - ^e 2 (l + cos0) 2 , (3.17) 



d cos 9 An 

where 9 is the angle between the ingoing v e and the outgoing z/ M , or equiva- 
lently, between the ingoing v e and the outgoing v^. Put another way, forward 
scattering is favored and backward scattering forbidden. This is due to angu- 
lar momentum conservation. The ingoing v e and v e have opposite helicities 
and, in the CM frame, opposite momenta, so that their combined spins add 
up to 1. The same is true for the outgoing particles so that backward scat- 
tering would violate angular momentum conservation. In the rest frame of 
the medium the ingoing u e tends to have energies of the order of its Fermi 
energy, while the ingoing v c tends to have energies of order T. Because for- 
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Fig. 3.7. — Differential and production rates d 2 n/dedt vs. neutrino 
energy for rj e — r] Ve — and T = 12 MeV. 
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ward scattering is favored, the outgoing tends to inherit the larger energy 
of the ingoing i/ e . 

The differences of the source spectra, however, do not translate into sig- 
nificant spectral differences of the and fluxes emitted from the SN 
core. While pair annihilations and nucleon bremsstrahlung are responsible 
for producing or absorbing neutrino pairs and thus their equilibration with 
the stellar medium below the "neutrino-energy sphere," other processes, no- 
tably u fl e ± scattering and nucleon recoils, are more efficient for the exchange 
of energy between neutrinos and the medium between the equilibration and 
transport spheres. In our numerical runs we will find in fact that adding the 
new process to a SN simulation primarily modifies the flux with only minor 
modifications of the spectrum. 

The numerical implementation is vastly simplified by exploiting the fact 
that z/ e and z/ e are in LTE in the region where the pair processes are effec- 
tive. Instead of using the actual neutrino distributions in the phase-space 
integration we can use the equilibrium distribution with the local medium 
temperature. This approximation breaks down at larger radii where the v e v e 
process is unimportant anyway. 

In order to further reduce computation time one of the remaining three 
phase-space integrations can be approximated by the analytic expressions 
given in Takahashi, El Eid, & Hillebrandt (1978). This also requires sim- 
plifying the blocking factors. With fi e = —fi e + > we can approximate 
the positron occupation number by a Maxwell-Boltzmann distribution. For 
He/T ^ 2 this holds to very good accuracy. The greatest deviation is at 
/x e /T = and yields blocking factors too low by about 10%. However, e~ 
and v e are always degenerate in the relevant regions. 



3.6 Scattering on Electrons and Electron 
Neutrinos 

Even though the scattering reactions on and on u e , u e are closely related 
to the pair annihilations discussed in the previous section, scattering on e ± 
is more important than on v e or z/ e . The matrix elements for these reactions 
are just the crossed versions of the leptonic pair processes, 

\M\ 2 = 8 G 2 F [(C V + C A fs 2 + (CV - C A fu 2 ] (3.18) 

spins 

with the weak interaction coefficients of Equations ()3.12j) or ()3.13j) for scat- 
tering on e~ or on z/ e , respectively. For s = 2k\ ■ ft2 and u = —2k\ ■ the 
momenta are assigned to the particles according to Figure 13.91 Crossing the 
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Fig. 3.9. — Leptonic scattering processes. 



matrix element Equation (j3.18j) again by interchanging u <-> t, we obtain 
scattering on e + or v e . This is also true for scattering of on e~ or z/ e ; 
scattering of on e + or u e brings us back to Equation (j3.18|) . The numerical 
procedure for calculating the rates is the same as for the pair-annihilation 
processes, given in the previous section. 

In contrast to the pair rates of the previous section, the differential scat- 
tering rates are monotonically rising functions of the degeneracy. In Fig- 
ure 13.101 we show the rates for scattering on e 1 * 1 (lower line) and the rates 
for scattering on z/ e or z/ e , normalized to the latter rate. Scattering on u e 
and z/ e is always more frequent than scattering on e in case of r) Uc = n e . How- 
ever, for realistic situations with r) Ve < r] e we expect that scattering on e 1 * 1 has 
1-2 times the rate of scattering on v e and v e . Therefore, neutrino-neutrino 
scattering is expected to be a relatively minor correction. In our numerical 
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Fig. 3.10. — Thermally averaged scattering rate for on e as a function 
of r/ e (lower line) and for on u e and u e (upper line) as a function of n Ue . 
The rates are normalized to the scattering rate on e 1 * 1 at r/ e = 0. We used 
T = 12 MeV and -q v = 0. 
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studies we will indeed find that this process has only a small effect on the 
neutrino spectra and fluxes. 

3.7 Other Reactions 

Interactions that also take place inside a SN, but are unimportant for our 
purpose, are absorption and emission of v e and z/ e by nuclei as well as scatter- 
ing of neutrinos on nuclei. In the vicinity of the neutrino sphere all nuclei are 
disintegrated. Reactions on nuclei take place at larger radii and are therefore 
important for the heating process behind the stalled shock wave. 

Throughout the past literature, oftentimes the plasmon process was taken 
into account as a source for neutrino pairs. The plasmon is an excitation of 
the medium that can emit a vv pair once it decays. In the relevant area 
close to the neutrino sphere this process is several orders of magnitude less 
frequent than the other pair processes. 



Chapter 4 



Setting the Stage for Neutrino 
Transport 

Proto-neutron star atmospheres are characterized by their radial profiles of 
temperature, nucleoli density, and the number of electrons per nucleon. In 
order to probe a variety of possible atmospheric structures we use two self- 
consistent profiles obtained by hydrodynamic simulations and, in addition, 
power-law parameterizations with constant electron fractions. After defining 
the optical depth and an effective mean free path for thermalizing processes 
we calculate the thermalization depths, i.e., the last points of interaction for 
all the thermalizing processes in each of the background models. This simple 
approach is a powerful test of the relative importance of various neutrino- 
matter interactions and their dependence on stellar parameters. 

4.1 Characterizing Proto-Neutron Stars 

All the neutrino interaction rates with the background medium depend on the 
density of scatterers and their temperature. Both density p and temperature 
T are functions of the radius. In the region relevant for neutrino spectra 
formation we assume that the proto-neutron star consist only of 

• protons 

• neutrons 

• electrons 

• positrons 

• neutrinos of all flavors. 
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In order to calculate every interaction rate for all radii, we need the number 
density and the chemical potential for each of these matter constituents as 
a function of radius in addition to the temperature. As mentioned earlier no 
muons and taus are present. Therefore, usually we set the chemical potential 
of to zero. Only in the case of weak magnetism a small degeneracy 
develops and we obtain a chemical potential for v^. 

Our neutron-star atmospheres are characterized by T(r), p(r), and the 
electron fraction per baryon Y c {r). Since m e <C m np , the density is given by 
the number densities of neutrons and protons times their masses 

p(r) = n n m n + n p m p . (4.1) 

In addition we use 

71 

Y e (r) = (4.2) 

Tin + U p 

to obtain the densities of neutrons and protons. By solving 

n n , P (r) = f de/ nj p(e,/x niP ) = f de — 6 M , (4.3) 

J J 1 + exp( — 

we obtain the chemical potentials of protons and neutrons /i n and /z p . 
Similarly, we determine p e by 

n B Y c = n e - -n e + = J de (/(e, /i e ) - /(e, -/i e )) . (4.4) 

A vanishing electron fraction does not mean that there are no electrons 
present, but rather that there is an equal thermal distribution of electrons 
and positrons. 

After all other chemical potentials are known, y, Ve is given by Equa- 
tion (I3~TH1) 

Pv e = Pe + - A^n • 

This completes the set of thermodynamic quantities needed for the transport 
of neutrinos in a SN core. 



4.2 Proto-Neutron-Star Profiles 

One possibility for our study of the influence of neutrino interactions and the 
stellar background model on the emerging neutrino spectra is to use proto- 
neutron star atmospheres that were obtained by self-consistent hydrodynamic 
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simulations. As a first example we use a model representative for the accre- 
tion phase; henceforth we will refer to it as the "Accretion-Phase Model I" 
(Fig. 14. lj) . It was provided to us by O. E. B. Messer and was already used in 
Raffelt (2001) for a more schematic study. Based on the Woosley & Weaver 
15 M Q progenitor model labeled sl5s7b, the Newtonian collapse simulation 
was performed with the SN code developed by Mezzacappa et al. (2001). 
The snapshot is taken at 324 ms after bounce when the shock is at about 
120 km, i.e., the proto- neutron star still accretes matter. In this simulation 
the traditional microphysics for transport was included, i.e., iso-energetic 
scattering on nucleons, e + e~ annihilation and i^e - scattering. 

As another self-consistent example (Accretion-Phase Model II, Fig. 14. 2j) 
we obtained a 150 ms post bounce model from M. Rampp (personal com- 
munication) that uses a very similar progenitor (sl5s7b2). The simulation 
includes an approximate general relativistic treatment in spherical symme- 
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Fig. 4.1. — Accretion-Phase Model I, a SN model 324 ms after bounce from 
a Newtonian calculation (O.E.B. Messer, personal communication). 
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try as described by Rampp & Janka (2002). The three neutrino flavors are 
transported with all relevant interactions except z/ e z/ e — ► v^v^ (Section 17.11 
and Rampp et al. 2002). 

When we use these self-consistent models with the same microphysics 
that was used in the original models by Messer and Rampp we find good 
agreement with their results. Changing the microphysical input then enables 
us to study the influence of the various interaction processes (Chapter EJ) on 
the emerging neutrino spectra in a systematic way, which has never been 
done before. Since we are simulating a static proto-neutron star and also 
vary the input physics, our simulations are not self consistent. However, once 
our findings are implemented in the self-consistent hydrodynamic simulations 
they qualitatively find the same results as we do (Buras et al. 2003b). 

The self-consistently obtained stellar background models are in the rele- 
vant range to a good approximation just power laws for the radial dependence 
of temperature and density. As a different approach we can thus parameterize 
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Fig. 4.2. — Accretion-Phase Model II, a SN core at 150 ms post bounce from 
a general-relativistic simulation. (M. Rampp, personal communication). 
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the profiles in a simple way and thereby study the influence of the stellar pro- 
files. Especially steep profiles that are characteristic for the late phase of the 
proto-neutron star are not available to us from self-consistent simulations. 
We use two power-law profiles of the form 

, = p„(^) p , r = T„(?)*, «. 5) 

with a constant electron fraction per baryon Y e . We adjust parameters such 
that (e) 20-25 MeV for the emerging neutrinos to obtain model atmo- 
spheres in the ball park of results from proto-neutron star evolution calcula- 
tions. We define a "steep" power-law model, corresponding to the one used 
by Raffelt (2001), and a "shallow" one; the characteristics are given in Ta- 
ble 4.1. The shallow model could be characteristic of a SN core during the 
accretion phase while the steep model is more characteristic for the neutron- 
star cooling phase. The constant electron fraction Y c is another parameter 
that allows us to investigate the relative importance of the leptonic processes 
as a function of the assumed Y c . 

Table 4.1. Characteristics of power-law models. 





Steep 


Shallow 
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q/p 
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p [10 14 gem" 3 ] 


2.0 
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T [MeV] 


31.66 


20.0 


r [km] 


10 


10 



4.3 Optical Depth vs. Thermalization Depth 

By comparing the optical depths due to different neutrino interactions we 
can assess the relative importance of the neutrino processes in the stellar 
background model. In principle the optical depth r is defined as the integral 
of the inverse mfp A -1 (see e.g. Shapiro & Teukolsky 1983, Suzuki 1989): 

roo | roo 

Ti(r,e)= dr \ = dr'a^e) m{r') , (4.6) 

Jr Ai[r , t) Jr 
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where % stands for any interaction channel, <7j is the corresponding cross 
section, and rii the number density. The optical depth specifies how frequently 
a specific process occurs on a neutrino's way out of the SN core. 

Since we want to know how efficient a certain process is at exchanging 
energy we have to take into account the qualitatively different nature of 
the neutrino-matter interactions. The relevant interaction processes fall into 
three classes: 

1. Creation processes, that are able to keep neutrinos in LTE. They 
can create and absorb neutrinos and therefore also contribute to the 
energy exchange between neutrinos and the background medium. 

2. Scattering with large energy transfer. In this case the neutrinos 
scatter on particles with a low mass compared to the neutrino's mo- 
mentum (i.e., ^ and neutrinos). 

3. Scattering on nucleons. Here the energy transfer is very small in a 
single scattering, The fact that allows for the effective mfp description 
discussed in the previous section. 

All three classes contribute to the opacity, but only the first two classes of 
reactions can exchange large amounts of energy in single interactions. 

Making use of this qualitative difference between nucleon scattering and 
the other interactions we can define the location (thermalization depth) 
where a neutrino of given energy last exchanged energy efficiently with the 
medium by a reaction such as v^er e~z/ M . Looking at the example of z/^e - 
scattering we have two different mfps, namely one for nucleon scattering At 
and one for the energy exchanging v^e~ process Xe- The total mfp is then 
given by A to t = l/(Ay 1 + A^ 1 ). Between two energy exchanging scatterings 
the neutrino undergoes many isoenergetic scatterings. This yields an effective 
mfp for the energy exchange reaction that is larger than \ E . 

The way of a neutrino through the medium can be considered a random 
walk. Therefore the traveled distance scales as the square root of time that is 
proportional to \^N, where N represents the number of interactions. For each 
reaction the probability of energy exchange is A to t/A.E, so the inverse is the 
number of reactions until the next energy exchange. Altogether the traveled 
distance from one energy-exchange reaction to the other corresponds to the 
effective mfp of the energy exchanging process: 




(4.7) 
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With the effective mfp we can then calculate the optical depth for thermal- 
ization 



When the optical depth becomes less than of order unity, the corresponding 
process is no longer relevant. The corresponding radius is referred to as the 
thermalization depth: 



where -Rtherm depends on the neutrino energy e. 

With this concept we can determine the energy-dependent sphere where 
a certain class of reactions becomes inefficient. For a specific neutrino energy 
the corresponding radius is given by the largest -Rtherm of any reaction in the 
class. For the first class this sphere corresponds to the number sphere, where 
particle creation freezes out, as defined in Section 12.51 The second class of 
reactions freezes out at the radius that determines the energy sphere. 

For the third class, i.e., neutrino- nucleon scattering, we are not able to 
obtain the thermalization depth. In each reaction the energy transfer is very 
small, thus many reactions are necessary to thermalize the neutrino. We can 
therefore not determine the location of the transport sphere, i.e., the radius 
at that neutrino- nucleon scattering becomes inefficient. 

As we show in the next section, due to the energy dependence it is not 
possible to define a unique neutrino sphere. Even for u e and i> e , where the 
/5-processes dominate and therefore the only relevant sphere is their num- 
ber sphere, its location is strongly energy dependent. The other flavors are 
emitted at their number sphere, lose energy in elastic scattering up to their 
transport sphere, and then diffuse before they start streaming freely from 
their transport sphere. Compared to v e and u e one would refer to the num- 
ber sphere as neutrino sphere of v^. The spectrum emitted at this neutrino 
sphere, however, is very different from the emitted spectrum. 

In order to illustrate the transport let us consider bremsstrahlung and 
iso-energetic scattering on nucleons only. For nucleon bremsstrahlung the 
energy dependence is rather weak and at the same time the density depen- 
dence is very strong. Therefore, to a good approximation one may picture 
the bremsstrahlung- number sphere as a blackbody surface that injects neutri- 
nos into the scattering atmosphere and absorbs those scattered back (Raffelt 
2001). With iso-energetic scattering on nucleons, the neutrino flux and spec- 
trum emerging from the transport sphere is then easily understood in terms 
of the energy- dependent transmission probability of the blackbody spectrum 




(4.8) 




2 

3' 



(4.9) 
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launched at the number sphere, that is in this simple picture equal to the 
energy sphere. The transport cross section scales as e 2 , implying that the 
transmitted flux spectrum is shifted to lower energies relative to the temper- 
ature at the energy sphere. This simple "filter effect" accounts surprisingly 
well for the emerging flux spectrum (Raffelt 2001). For typical conditions 
the mean flux energies are 50-60% of those corresponding to the blackbody 
conditions at the energy sphere. 

Moreover, it is straightforward to understand that the effective tempera- 
ture of the emerging flux spectrum is not overly sensitive to the exact location 
of the number sphere. If the pair processes are somewhat more effective, the 
number sphere is at a larger radius with a lower medium temperature. How- 
ever, the scattering atmosphere has a smaller optical depth so that the higher- 
energy neutrinos are less suppressed by the filter effect, partly compensating 
the smaller energy-sphere temperature. For typical situations Raffelt (2001) 
found that changing the bremsstrahlung rate by a factor of 3 would change 
the emerging neutrino energies only by some 10%. This finding suggests that 
the emitted average neutrino energy is not overly sensitive to the details of 
the pair processes. 

With the complete set of interactions neutrinos undergo energy exchang- 
ing scattering reactions below the energy sphere before they enter the dif- 
fusive regime. Therefore the mean energy of the emerging flux is even less 
sensitive to the exact location of the number sphere than found by Raffelt 
(2001). In Section IB~3*1 we show that changing the bremsstrahlung rate by a 
factor of three up or down has almost no effect on the mean energy of the 
emerging neutrino flux. 

4.4 Thermalization Depths in Our Stellar 
Models 

We now turn to applying the concept of thermalization depths to the neutron- 
star atmospheres described above (Section l4.2j) . We consider the neutrino mfp 
for nucleon bremsstrahlung NN — > NNV^z/^, pair annihilation e + e~ — > z^P^ 
and v e v e — > v^v^ and scattering on charged leptons — > e ± u fl . The 
reaction rates used were described in Chapter El 

In Figures 14.31 and 14.41 we give the thermalization depth -Rtherm as a 
function of neutrino energy e for the two hydrodynamically self-consistent 
accretion-phase models. From top to bottom the panels show the results 
for u e , i> e , and respectively. The step-like curves represent the tempera- 
ture profiles in terms of the mean neutrino energy, (e„) = 3.15 T for non- 
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degenerate neutrinos at the local medium temperature; the steps correspond 
to the radial zones of our numerical simulation. Slightly above these steps the 
thin line represents the mean energy of neutrinos flowing in radial direction 
obtained by our numerical simulation. The other curves represent -Rtherm for 
bremsstrahlung (b), e + e~ annihilation (p), u c u e annihilation (n), and scat- 
tering on e (s). In the case of u e and u e we do not include bremsstrahlung 
and v e v e annihilation. Particle creation is dominated by the charged current 
reactions on nucleons (urea). 

For the power-law models we show -Rtherm for in Figures 14.51 and 14.61 
The different panels correspond to the indicated values of the electron fraction 
Y e . Note that Y e represents the net electron density per baryon, i.e., the e~ 
density minus that of e + so that Y e = implies that there is an equal thermal 
population of e~ and e + . 

The absorption rate for the bremsstrahlung process varies approx- 
imately as e _1 , the z/^N transport cross section as e 2 so that the inverse 
mfp for thermalization varies only as e 1 / 2 . This explains why -Rtherm f° r 
bremsstrahlung is indeed quite independent of e. Therefore, bremsstrahlung 
alone allows one to specify a rather well-defined number sphere. The other 
processes depend much more sensitively on e so that a mean energy sphere 
is much less well defined. 

Both electron scattering and the leptonic pair processes are so ineffective 
at low energies that true LTE can not be established even for astonishingly 
deep locations. Bremsstrahlung easily "plugs" this low-energy hole so that 
one can indeed expect LTE for all relevant neutrino energies below a cer- 
tain radius. For higher energies, the leptonic processes dominate and shift 
the energy sphere to larger radii than bremsstrahlung alone. The relative 
importance of the various processes depends on the density and temperature 
profiles as well as Y e . 

To assess the role of the various processes for the overall spectra forma- 
tion one needs to specify some typical neutrino energy. One possibility would 
be (e) for neutrinos in LTE. Another possibility is the mean energy of the 
neutrino flux, in particular the mean energy of those neutrinos which actu- 
ally leave the star. In Figures f4.3H4. 61 the thin line represents this flux mean 
energy. In the case of it is clearly visible that even after the last energy 
exchanging process, i.e., e scattering, became ineffective the mean energy 
of the flux still drops. This drop is due to scattering on nucleons. Among the 
pair processes it is always the u e P e annihilation process that has the largest 
-R t h er m in the range of flux energy. However, especially in the steep power-law 
profile the intersection of the bremsstrahlung' s and v e v e annihilation's ther- 
malization depths are rather close at the mean flux energy. The e + e~ process 
is always less important than u e P e , as was already discussed in Section 13*31 
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In our numerical simulations we find the same results (Chapter EJ) . 

The final shape of the z/ M spectra and therefore also their mean energy 
cannot depend much on any of the ^-creation processes. The energy sphere 
represented by -Rtherm of e ± scattering stretches in all profiles far beyond the 
other thermalization depths. Additionally, also the scattering on nucleons 
has a significant impact on the mean energies. Its importance compared to 
electron scattering is very difficult to guess. Recall, it is not possible to define 
a thermalization depth for nucleon recoils. Both processes are qualitatively 
very different in that neutrinos transfer only a small fraction of their energy 
per nucleon scattering. Numerically we have compared the importance of all 
processes and present our results in Chapter^ 
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R therm t km ] 

Fig. 4.3. — -Rtherm as a function of neutrino energy e for our Accretion-Phase 
Model I. From top to bottom the panels show the results for u e , u e , and 
v^. Energy exchanging processes: bremsstrahlung (solid line), e + e~ annihila- 
tion (dashed), u e u e annihilation (dotted), and scattering on e ± (dash-dotted). 
"Urea" denotes the charged- current reaction of u e and u e on nucleons. The 
steps represent (e) = 3.15 T; the thin line close to 3.15 T is the mean energy 
of the flux (e)fl ux , cf. Chapter El 
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Fig. 4.4. — Same as Fig 



;ure 14.31 for the Accretion-Phase Model II. 
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Fig. 4.5. — -Rtherm for ^ in the steep power-law model with the indicated 
values of Y e . This figure corresponds to the bottom panel of Figure 1431 
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Fig. 4.6. — Same as Fig. 14.51 for the shallow power-law model. 



Chapter 5 

Characterizing Non-Thermal 
Neutrino Spectra 

The neutrino fluxes emerging from a SN core are not thermal. There are nu- 
merous ways to describe spectral characteristics with a few parameters, for 
example energy moments. For many applications it is useful to have an ana- 
lytic description of the spectra. Throughout the literature it is common to use 
an nominal Fermi-Dirac distribution. We present a simpler two-parameter fit 
that describes our Monte Carlo data better, and compare both descriptions. 



5.1 Global Parameters 

While a SN core can be crudely understood as a blackbody source for neu- 
trinos of all flavors, the escaping fluxes are not strictly thermal. A number 
of observables have been defined in the literature in order to describe the 
properties of these spectra. 

Let us assume f(e,/i) is the distribution of neutrinos as a function of 
energy e and /i = cos 9, where 9 is the angle between the radial direction and 
the direction of motion. Far away from the SN core neutrinos stream freely 
in radial direction, therefore — 

As the most intuitive quantity characterizing the spectrum the mean en- 
ergy is commonly used, 

^ = I£dee /+ 1 dfi fje, fx) 
J^de f+l dfi /(e, fi) 

In this ratio the numerator is the energy density and the denominator is the 
number density. 
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,\ Jo u " J-i u ^FJ^.w fK > 

( e )fl^ = roc ~ F+i ~ .- ^t: ~ • ( 5 - 2 ) 



Within the number sphere neutrinos are in LTE, but in the outer regions 
some net flux of neutrinos develops. Therefore, it is often useful to extract 
spectral characteristics for those neutrinos which are actually flowing by re- 
moving the isotropic part of the distribution. Specifically, we define the mean 
flux energy by 

f^dee J^dfifj,f(e,fjL) 

Outside the transport sphere all neutrinos will flow essentially in the radial 
direction, implying that the angular distribution becomes a delta-function in 
the forward direction so that (e)fl ux = (e). In the trapping regions the two 
averages are very different because the distribution function is dominated by 
its isotropic term. However, for most discussions only the emerging spectra 
are relevant and both definitions of the mean energy, Equations (jQJ) and 
(15. 2 j) are equivalent. The right hand side of Equation (|5.2|) then corresponds 
to the ratio of luminosity to number flux. 

Deep inside the core, below the number sphere neutrinos are in LTE. 
With a vanishing chemical potential their distribution function is given by 

f{e^)= e \ (5.3) 
1 + exp(e/T) 

and therefore 

<e) = ^r«3.1514T. (5.4) 

We use this relation to compare the numerically obtained mean energy to 
the medium temperature, cf. e.g. Figure PPl 

Since the emerging spectra are not thermal, additional characteristics 
are useful. A rather general way to characterize the spectrum beyond mean 
values is by its energy moments (Janka & Hillebrandt 1989a) 



f Q °°def+tdixf(e,ri 



<o = J0 r r;. X 7. • ( 5 - 5 ) 



Usually, values are given only up to n=2 in the literature, i.e., the mean 
energy and the second moment. 

With the first two moments an easy way to see by how much a spectrum 
deviates from being thermal is the pinching parameter, defined by Raffelt 
(2001). It is basically the ratio of the first two moments 



(5.6) 
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normalized such that a Fermi-Dirac distribution with vanishing chemical po- 
tential yields p = 1. Therefore, 

(e 2 ) FD 486000 C3C5 , v 

For a Maxwell-Boltzmann distribution this quantity would be 4/3. A spec- 
trum that is thermal up to its second moment has p = 1, while p < 1 signi- 
fies a pinched spectrum (high-energy tail suppressed), p > 1 an anti-pinched 
spectrum (high-energy tail enhanced). 

In some publications the root-mean-square energy (e) rms is given instead 
of the average energy. The rms-energy involves the ratio of third and first 
energy moment 



I °°def+tdu.e3f(e,n) 



\ Jo°def±}dtief{e,fi) \ (e) 



" (5. 



This characteristic spectral energy is useful for estimating the energy transfer 
from neutrinos to the stellar medium in reactions with cross sections propor- 
tional to e 2 . However, comparing the mean energy defined in Equation ()5.1|) 
and the rms-energy is difficult, because the relation of both quantities de- 
pends on the spectral shape. 



5.2 Analytic Fits 

For applications like neutrino oscillations it is useful to have an analytic fit to 
the SN neutrino spectra. We present a two-parameter fit that reproduces the 
first five or more moments of our numerically obtained spectra to a very good 
accuracy (Section I7.3J) . In addition to the overall normalization we introduce 
the parameters a and e by defining 

/ Q (e) = i(0V (Q+1)e/ % (5.9) 

with the normalization 

c = de/ a (e) = (a + l)- (a+1) r(a + l)e. (5.10) 
Jo 

The average energy is always (e) = e, while a controls the pinching of the 
function. Just like the mean energy all observables discussed in the previous 
section can be expressed in terms of a and e by simple analytic relations. 
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Fig. 5.1. — The ratio of rms energy and mean energy. The left panel shows 
the dependence on a as given in Equation (|5.14j) ; in the right panel we plot 
the dependence on r\. 



The pinching parameter p is closely related to the width 



w = s/W) - (e) 2 = y/ap-1 (e) = e , (5.11) 

with a ~ 1.3029 as defined in Equation ()5.7j) . The relation for the pinching 
parameter can be read off Equation ()5.11|) . 

{e) z 1 + a 

For the more general case of the ratios of arbitrary neighboring moments 
we find 

= (5-i3) 

(e 71 " 1 ) 1 + a 

that can be easily generalized to arbitrary ratios of moments with Equa- 
tion flnnD. 

Taking the ratio of third and first energy moment yields the rms-energy 



(4 



As shown in in the left panel of Figure 15.11 (e) rms is always some 25-45% 
larger than (e) in the range of a = 2-5, that turns out to be the parameter 
range obtained by numerical simulations. Therefore, the two different energy 
averages contain information on the spectral shape. 
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From any three moments it is possible to obtain the normalization, a, 
and e. From all our studies we found that the spectra are very well described 
by these three parameters. 

In the literature one frequently encounters a nominal Fermi-Dirac dis- 
tribution characterized by a temperature T and a degeneracy parameter r\ 
according to 

m = ^ V (5-15) 



1 + exp (~ 



n 



as a representation of the SN neutrino spectrum ( Janka & Hillebrandt 1989a) . 
The functional form is motivated by the equilibrium distribution of neutrinos 
inside the star. The values of the parameters T and rj, however, are not easy 
to interpret, because the emerging spectra are not thermal. 

The functional form of f v is more complicated than f a . Most calculations 
involving f n have to be carried out numerically and/or approximately. In 
Figure 15.21 we show (e/T) and p as a function of rj. Up to second order, 




Fig. 5.2. — Mean energy and pinching parameter as a function of the degen- 
eracy parameter for a Fermi-Dirac distribution. As dashed lines we show the 
expansions given in Equation (|5.16jl . 
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expansions are 

(e/T) « 3.1514 + 0.1250 n + 0.0429 rj 1 

p « 1 - 0.0174//- 0.0046 r/ 2 . (5.16) 

These expansions are shown in Figure IB~2l as dashed lines. 

There are no analytic expressions for the energy moments and therefore 
also the rms-energy can only be obtained numerically. For the special case of 
rj = the expression is 



^ rms = V 441 71 T ~ 4,5622 T ' (5 ' 1?) 

With Equation f)5.4|) this corresponds to (e) « 0.691 (e) rms , but recall that 
SN neutrino spectra are not thermal spectra and rj ^ 0. The numerically 
obtained dependence of the ratio of rms and mean energy is displayed in the 
right panel of Figure l5~Tl It is rather similar to the case of the a fit. 

In view of analytic simplicity the a fit is certainly superior, but the shapes 
of both functions are rather similar. The numerically determined spectra are 
similarly well described by both f v and f a . As we show in the following for 
a reasonable range for the parameters both functions differ by no more than 
about 10%. 

In the upper panel of Figure I5~3l we show f a (e), the integral normalized 
to unity, for several values of a. The broadest curve is for a = 2 while for 
the narrower ones the width given in Equation ()5.11j) was decreased in 10% 
decrements as shown in Table 15.11 The middle panel of Figure 15.31 shows 
the corresponding curves /^(e) with the n- values given in Table 15.11 The 
broadest curves in each panel are identical and correspond to e 2 exp(— 3e/e) 
with a width w = (e) / \/3. 

Table 5.1. Parameters for fit-functions of Figure l5~?fl 



Width 


a 


V 


wo = (e)/V3 


2. 


— oo 


0.9 w 


2.7037 


1.1340 


0.8 w 


3.6875 


2.7054 


0.7w 


5.1225 


4.4014 


0.6 w 


7.3333 


6.9691 


0.5 Wq 


11. 


13.892 
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Fig. 5.3. — Normalized fit functions. Upper panel: a-fit according to Equa- 
tion (|5.9|) . Middle panel: Fermi-Dirac fit according to Equation (|5.15|) . In 
both panels the broadest curve corresponds to /(e) = e 2 exp(— 3e/e), i.e., to 
a = 2 and rj = — oo, respectively. For the other curves the width was de- 
creased in decrements of 10%, see Table l5~Tl Bottom panel: Ratio of the fits 
fa/ ft] f° r the first four cases. 



56 



Chapter 5. Characterizing Non-Thermal Neutrino Spectra 



The limiting behavior of / a (e) for large a is S(e — e) while for f n {e) the 
limiting width is wo/\/h ~ 0.44721 wq. Evidently the curves f a (e) can ac- 
commodate a much broader range of widths than the curves /^(e). 

We find that the neutrino spectra are always fit with parameters in the 
range 2<a<5or0<7/<4 (Section l7.1j) . i.e. with a width above about 
0.75 wq. In the bottom panel of Figure 15.31 we show the ratios of the fit 
functions for the widths down to 0.7 wq. From this plot we find that except for 
the lowest energies and very high energies the two fit functions are equivalent 
to better than 10%. Therefore, the two types of fits are largely equivalent for 
most practical purposes. 

On the basis of a few high-statistics Monte Carlo runs we show in Sec- 
tion 17.31 that the numerical spectra are actually better approximated over 
a broader range of energies by the "power-law" fit functions / a (e). In addi- 
tion, these functions are more flexible at representing the high-energy tail of 
the spectrum that is most relevant for studying the Earth effect in neutrino 
oscillations. 



Chapter 6 



Relative Importance of 
Neutral-Current Reactions 

We briefly introduce our Monte Carlo code, that was developed to study 
neutrino transport in proto-neutron-star atmospheres. For the proto-neutron- 
star atmospheres presented in Section l4~2l we show our results that were ob- 
tained by switching on and off various neutral-current interaction processes. 

The traditional set of processes included in transport turns out to be 
a rather poor approximation. Contrary to what had been assumed in the 
past, nucleon-nucleon bremsstrahlung and the z/ e z/ e pair process are the main 
sources for v^p^ production. Above the number sphere the spectra are shaped 
by electron scattering and scattering on recoiling nucleons. The traditional 
e + e~ — > v e v e process and scattering on v e jv e turn out to be rather negligible. 

These results and those of Chapter represent a major part of this dis- 
sertation and were already published in a similar form as: 
M. T. Keil, G. G. Raffelt, and H. T. Janka, "Monte Carlo study of supernova 
neutrino spectra formation," Astrophys. J. in press, astro-ph/0208035, 

6.1 Monte Carlo Setup 

For the stellar background models introduced in Section 14.21 we have run 
our Monte Carlo code using various sets of neutrino-matter interactions. In 
Appendix |B] we give a detailed description of the code. Our main interest 
is to assess the impact of the various neutrino interactions on the flux and 
spectrum formation. Since neutrinos are in LTE below the number sphere, 
it is sufficient to simulate the neutrino transport starting slightly below that 
radius where we have to specify a boundary condition. 

At that lower boundary of our proto-neutron-star atmosphere we always 
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impose a blackbody boundary condition. We assume neutrinos to be in LTE, 
i.e., at the local temperature and the appropriate chemical potential. For nu- 
merical runs discussed in this chapter the latter is taken to vanish, because 
only are involved and weak magnetism is discussed in the next chapter. As 
a consequence of this inner boundary condition, the luminosity emerging at 
the surface is generated within the computational domain and calculated by 
our Monte Carlo transport. A small flux across the inner boundary develops 
because of the negative gradients of temperature and density in the atmo- 
sphere. Its magnitude depends on the radial resolution of the proto-neutron- 
star atmosphere and is therefore an artifact of the numerical implementation. 
However, for our setups this flux is small compared to the luminosity at the 
surface and therefore the emerging neutrino spectra do not depend on the 
lower boundary condition. 

The shallow energy dependence of the thermalization depth of the nu- 
cleon bremsstrahlung implies that whenever we include this process it is not 
difficult to choose a reasonable location for the lower boundary. Taking the 
latter too deep in the star is very CPU-expensive as one spends most of the 
simulation for calculating frequent scatterings of neutrinos that are essen- 
tially in LTE. 

The outer boundary is determined by the requirement that neutrinos 
stream freely. In the outer regions the transport is not CPU-expensive. There- 
fore, the exact location of the outer boundary is not very crucial. However, 
we always use a grid of 30 equally spaced radial zones and can thus obtain 
better spatial resolution by choosing the outer boundary fairly low. 

We always include u^N scattering as the main opacity source. For energy 
exchange, we switch on or off the processes given in Table 16.11 We never 
include inelastic nucleon scattering ^NN — > NNz/ M as this process is never 

Table 6.1. Neutral-current processes included in our tranport. 



Process 


Label 


bremsstrahlung 


b 


nucleon recoil 


r 


scattering on e + and e~ 


s 


e + e~ pair annihilation 


P 


v e v e annihilation 


n 


scattering on u e and v e 


sn 
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important relative to recoil (Raffelt 2001). For confirming that scattering on 
u e and u e is really unimportant if is included (Section l3.6|) we used the 
Accretion-Phase Model I. In all other models we ignore this process. We also 
neglect v^v^ or v^v^ scattering even though such processes may have a larger 
rate than some of the included leptonic processes. Processes of this type do 
not exchange energy between the neutrinos and the background medium. 
They are therefore not expected to affect the emerging fluxes and should 
also have a minor effect on the emitted spectra. We study weak magnetism 
in Chapter 

6.2 Accretion-Phase Model I 

Our first goal is to assess the relative importance of different energy-exchange 
processes for the transport. As a first example we begin with our Accretion- 
Phase Model I. The results from our numerical runs are summarized in Ta- 
ble IHI21 where for each run we give (e)fl ux , (e 2 )fl ux , our fit parameter a deter- 
mined by Equation (|5.12j) . and the pinching parameter pfl ux for the emerging 
flux spectrum, defined in Equation (|5.fij) . We fitted the temperature and 
degeneracy parameter of an effective Fermi-Dirac spectrum producing the 
same first two energy moments, and give the luminosity. Only (e) flux , (e 2 )flux, 
and the luminosity were obtained from the numerical simulations. All other 
characteristics were calculated from these first two energy moments. 

The first row contains the muon neutrino flux characteristics of the orig- 
inal Boltzmann transport calculation by Messer. To make a connection to 
these results we ran our code with the same input physics, i.e., scatter- 
ing (s) and e + e~ annihilation (p). There remain small differences between the 
original spectral characteristics and ours. These can be caused by differences 
in the implementation of the neutrino processes, by the limited number of 
energy and angular bins in the Boltzmann solver, the coarser resolution of 
the radial grid in our Monte Carlo runs, and by our simple blackbody lower 
boundary condition. We interpret the first two rows of Table 16.21 as agreeing 
sufficiently well with each other that a detailed understanding of the differ- 
ences is not warranted. Henceforth we will only discuss differential effects 
within our own implementation. 

In the next row (bsp) we include nucleon bremsstrahlung which has the 
effect of increasing the luminosity by a sizable amount without affecting 
much the spectral shape. This suggests that bremsstrahlung is important as 
a source for v^v^ pairs, but that the spectrum is then shaped by the energy- 
exchange in scattering with e^. In the next row (bp) we switch off e ± scatter- 
ing so that no energy is exchanged except by pair-producing processes. The 
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spectral energy indeed increases significantly. However, the biggest energy- 
exchange effect in the scattering regime is nucleon recoil. In the next two 
rows we include recoil (brp) and then additionally scattering (brsp), both 
lowering the spectral energies and also the luminosities. 

The picture of all relevant processes is completed by adding v e v e pair 
annihilation (brspn), which is similar to e + e~ pair annihilation, but a factor 
of 2-3 more important ( Section I3.5|) . The luminosity is again increased, an 
effect which is understood in terms of our blackbody picture for the number 
and energy spheres. In the lower panel of Figure 14.31 we see that -Rtherm 
moves to larger radii once "n" is switched on, the radiating surface of the 
"blackbody" increases and more pairs are emitted. For both "p" and "n" 
-Rtherm is strongly energy dependent and therefore it is impossible to define a 
sharp thermalization radius. Switching off "r" again (bspn) shows that also 
with "n" included, "r" really dominates the mean energy and shaping of the 
spectrum. 



Table 6.2. Monte Carlo results for Accretion-Phase Model I. 



Enerj 


s;y exchan. 






(e 2 )flux 


a 


Pflux 


T 


T) 


K 


original run 




17.5 


388. 


2.7 


0.97 


5.2 


1.1 


14.4 




s 


P 




16.6 


362. 


2.2 


1.01 


5.3 


-0.3 


15.8 


b 


s 


P 




16.3 


351. 


2.1 


1.02 


5.4 


-2.2 


19.1 


b 




P 




17.8 


419. 


2.1 


1.02 


5.9 


-1.9 


20.1 


b r 




P 




15.1 


285. 


3.0 


0.96 


4.3 


1.6 


18.6 


b r 


s 


P 




14.2 


255. 


2.8 


0.98 


4.2 


1.1 


14.8 


b r 


s 


P 


11 


14.4 


264. 


2.7 


0.97 


4.3 


1.2 


17.6 


b 


s 


P 


11 


16.6 


358. 


2.3 


1.00 


5.2 


0.2 


21.7 




s 


P 


11 


16.9 


369. 


2.4 


0.99 


5.3 


0.4 


20.2 


b r 


s 






14.0 


251. 


2.6 


0.99 


4.3 


0.6 


13.1 


b r 


s 




11 


14.4 


263. 


2.7 


0.97 


4.3 


1.2 


17.0 


- r 


s 


P 




14.5 


265. 


2.8 


0.97 


4.3 


1.2 


13.0 


- r 


s 


P 


11 


14.7 


269. 


3.1 


0.96 


4.2 


1.7 


16.8 


b r 


sn 


P 


11 


14.3 


260. 


2.7 


0.97 


4.3 


1.2 


17.9 



Note. — The Energy exchange is specified in Table 16.11 We give 
(e)flux and T in MeV, (e )fl ux in MeV 2 , and L u in 10 51 erg s -1 . 
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To study the relative importance of the different pair processes, we switch 
off the leptonic ones (row "brs") and compare this to only the leptonic pro- 
cesses (row "rspn"). In this stellar model both types contribute significantly. 
Comparing then "brsp" with "brsn" shows that among the leptonic processes 
"n" is clearly more important than "p" . 

The last row (brsnpn) includes in addition to all other processes scattering 
on v e and u c . As already shown in Section |33J this process is about half as 
important as scattering on and its influence on the neutrino flux and 
spectra is negligible. We show this case for completeness but do not include 
scattering on v e and i> e for any of our further models. 

In order to illustrate some of the cases of Table 16.21 we show in the upper 
panel of Figure 16.11 several flux spectra from high-statistics Monte Carlo 
runs. Starting again with the input physics of the original hydrodynamic 




Fig. 6.1. — High-statistics spectra for Accretion-Phase Model I with different 
input physics as in Table Upper panel: Differential particle fluxes. Lower 
panel: Spectra normalized to equal particle fluxes. 
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simulation (sp) we add bremsstrahlung (b), recoil (r), and finally z/ e z/ e pair 
annihilation (n). Each of these processes has a significant and clearly visible 
influence on the curves. The pair-creation processes ("b" and "n") hardly 
change the spectral shape but increase the number flux, whereas recoil (r) 
strongly modifies the spectral shape. In the lower panel of Figure lfi~D we show 
the same curves, normalized to equal particle fluxes. In this representation 
it is particularly obvious that the pair processes do not affect the spectral 
shape. 

The very different impact of pair processes and nucleon recoils has a 
simple explanation. The thermalization depth for the pair processes is deeper 
than that of the energy-exchanging reactions, i.e., the number sphere is below 
the energy sphere. Therefore, the particle flux is fixed more deeply in the 
star while the spectra are still modified by energy-exchanging reactions in 
the scattering atmosphere. 

6.3 Steep Power Law 

As another example we study the steep power-law model defined in Equa- 
tion (|4.5jl and Table 4.1. This model is supposed to represent the outer 
layers of a late-time proto-neutron star but without being hydrostatically 
self-consistent. It connects directly with Raffelt (2001), where the same pro- 
file was used in a plane parallel setup, studying bremsstrahlung and nucleon 
recoil. The results of our runs are displayed in Table 16.31 and agree very 
nicely with those obtained by Raffelt (2001), corresponding to our cases "b" 
and "br". 

For investigating the importance of leptonic processes, we run our code 
with a variety of neutrino interactions and in addition assume a constant 
electron fraction Y e throughout the whole stellar atmosphere. This assump- 
tion is somewhat artificial, but gives us the opportunity to study extreme 
cases in a controlled way. In the relevant region Y e = 0.5 yields the highest 
possible electron density. In addition we study the electron fraction being 
one order of magnitude smaller, Y e = 0.05, and finally the extreme case with 
an equal number of electrons and positrons, Y e — 0. 

The first leptonic process we consider is e + e~ pair annihilation. Compar- 
ing the rows "bp" with the row "b" shows a negligible effect on the spectrum, 
but a rise in luminosity. Increasing Y e brings the luminosity almost back to 
the "b" case, because the electron degeneracy rises and the positron density 
decreases so that the pair process becomes less important. 

Adding scattering on forces the transported neutrinos to stay closer 
to the medium temperature, i.e., reduces their mean energy. Of course, the 



6.3. Steep Power Law 



63 



scattering rate increases with the number density of electrons and positrons, 
i.e., for higher Y e we get lower spectral energies. For the luminosity the sit- 
uation is more complicated. Since the neutrino flux energies decrease when 
we switch on scattering we would expect a lower luminosity. However, the 
opacity of the medium to neutrinos is strongly energy dependent and low 
energy neutrinos can escape more easily than high-energy ones, increasing 
the number flux. On balance, the "bsp" luminosities are larger compared to 
the "bp" ones. 

To compare the scattering on e 1 * 1 with that on nucleons, we turn off "s" 
again and instead switch on recoil (r). Qualitatively, the energy exchange is 
very different from the earlier case. In the scattering on e 1 * 1 a neutrino can 
exchange a large amount of energy, while for scattering on nucleons the energy 
exchange is small. But since neutrino-nucleon scattering is the dominant 
source of opacity that keeps the neutrinos inside the star, the scatterings are 
very frequent. This leads to a stronger suppression in the high-energy tail of 
the neutrino spectrum and therefore to a visibly smaller mean flux energy 



Table 6.3. Monte Carlo results for the steep power-law model. 



Energy exchange 




Y c 


(e)flux 


(e 2 )flux 


a 


Pflux 


T 


V 




b 












25.8 


962. 


1.2 


1.11 






21.0 


b 


r 










19.5 


487. 


2.6 


0.98 


6.0 


0.7 


14.5 


b 






p 







25.4 


890. 


1.6 


1.06 






23.8 


b 






p 




0.05 


25.6 


908. 


1.6 


1.06 






23.2 


b 






p 




0.5 


25.5 


917. 


1.4 


1.08 






21.6 


b 




s 


p 







24.2 


787. 


1.9 


1.03 






24.5 


b 




s 


p 




0.05 


23.8 


753. 


2.0 


1.02 






24.5 


b 




s 


p 




0.5 


21.3 


591. 


2.3 


1.00 


6.8 


-0.3 


23.1 


b 


r 




p 




0.5 


20.0 


507. 


2.7 


0.98 


6.0 


1.0 


16.8 


b 


r 


s 


p 







20.3 


518. 


2.9 


0.97 


5.9 


1.4 


19.7 


b 


r 


s 


p 




0.05 


20.3 


518. 


2.9 


0.97 


5.9 


1.4 


19.5 


b 


r 


s 


p 




0.5 


19.6 


488. 


2.7 


0.98 


5.9 


1.1 


18.7 


b 


r 


s 


p 


n 





20.7 


535. 


3.0 


0.96 


5.8 


1.8 


23.9 


bx3 


r 


s 


p 


n 


0.05 


20.3 


522. 


2.7 


0.97 


6.0 


1.3 


24.2 


b 


r 


s 


p 


n 


0.05 


20.6 


530. 


3.0 


0.96 


5.9 


1.7 


23.8 


bx0.3 


r 


s 


p 


n 


0.05 


20.7 


534. 


3.1 


0.96 


5.8 


1.8 


23.4 


b 


r 


s 


p 


n 


0.5 


19.8 


499. 


2.7 


0.97 


5.9 


1.2 


21.4 
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and lower effective spectral temperature, but higher a and higher effective 
degeneracy. Many nucleon scatterings, however, are needed to downgrade the 
high-energy neutrinos (different from scattering). Therefore neutrinos stay 
longer at high energies and experience a larger opacity and a larger amount 
of backscattering. This suppresses the neutrino flux significantly. 

In the runs including both scattering reactions (brsp), we find a mixture 
of the effects of e 1 * 1 and nucleon scatterings and an enhanced reduction of the 
mean flux energy. 

Finally, adding the neutrino pair process yields almost no change in en- 
ergy and pinching, but an increased luminosity as expected from the analo- 
gous case in Section l6~21 Although this profile is rather steep, leptonic pair 
processes are still important (Figure PO>j) . 

In order to estimate the sensitivity to the exact treatment of nucleon 
bremsstrahlung we have performed one run with the bremsstrahlung rate 
artificially enhanced by a factor of 3, and one where it was decreased by a 
factor 0.3. All other processes were included. The emerging fluxes and spectra 
indeed do not depend sensitively on the exact strength of bremsstrahlung as 
argued in Section l4~3l 

6.4 Shallow Power Law 

For the shallow power law (Table 16 .4|) almost the same discussion as for the 
steep case applies. As we can already infer from Figure WJ)\ leptonic processes 
are more important. This leads to a much higher increase of the neutrino flux 
once "p" or "n" are included, and to stronger spectral pinching when e + e~ 
annihilation is switched on. Scattering on downgrades the transported 
neutrino flux by a larger amount. 

6.5 Summary 

We find that the spectra are reasonably well described by the simple pic- 
ture of a blackbody sphere determined by the thermalization depth of the 
nucleonic bremsstrahlung process, the "filter effect" of the scattering atmo- 
sphere, and energy transfers by nucleon recoils. This is also true for the 
flux in case of steep neutron-star atmospheres. For more shallow atmospheres 
pair annihilation (e + e~ and v e Ve), however, yields a large contribution to the 
emitted flux and e 1 * 1 scattering reduces the mean flux energy significantly. 
It is therefore important for state-of-the art transport calculations to include 
these leptonic processes. The traditional process e + e~ — > v^p^ is always sub- 
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dominant compared to v c v c — * v^v^. The relative importance of the various 
reactions depends on the stellar profile. 

Neutrinos emitted from a blackbody surface and filtered by a scatter- 
ing atmosphere without recoils and leptonic processes have an anti-pinched 
spectrum (Raffelt 2001). However, after all energy-exchanging reactions have 
been included we find that the spectra are always pinched. When described 
by our a-fit, a ranges from about 2.5 to 3.5. For an effective Fermi-Dirac 
distribution, the nominal degeneracy parameter r\ is typically in the range 
1-2, depending on the profile and electron concentration. 

The traditional set of neutrino interactions that is commonly included 
in hydrodynamic SN simulations is not sufficient to account for spec- 
tra formation. Nucleon recoils significantly lower the mean energy, and 
bremsstrahlung as well as the v e v c pair process cause higher luminosities 
compared to results obtained with the traditional interactions. 



Table 6.4. Monte Carlo results for the shallow power-law model. 



Energy exchange 




(e)flux 


(e 2 )flux 


a 


Pflux 


T 


7] 


L v 


b 










27.7 


1120. 


1.2 


1.12 






20.3 


b r 










20.1 


521. 


2.5 


0.99 


6.3 


0.4 


13.4 


b 




P 







27.7 


974. 


2.7 


0.98 


8.3 


1.0 


43.1 


b 




P 




0.05 


27.9 


990. 


2.7 


0.98 


8.3 


1.1 


43.3 


b 




P 




0.5 


28.3 


1019. 


2.7 


0.98 


8.5 


1.0 


38.3 


b 


s 


P 







25.5 


830. 


2.6 


0.98 


7.6 


1.1 


46.2 


b 


s 


P 




0.05 


25.4 


815. 


2.8 


0.97 


7.5 


1.2 


46.3 


b 


s 


P 




0.5 


23.5 


706. 


2.6 


0.98 


7.1 


1.0 


44.8 


b r 




P 




0.5 


22.5 


624. 


3.3 


0.95 


6.1 


2.2 


33.1 


b r 


s 


P 







22.3 


612. 


3.3 


0.95 


6.1 


2.1 


39.6 


b r 


s 


P 




0.05 


22.2 


609. 


3.2 


0.95 


6.1 


2.1 


39.1 


b r 


s 


P 




0.5 


21.7 


585. 


3.1 


0.95 


6.1 


1.9 


39.2 


b r 


s 


P 


n 





22.2 


608. 


3.3 


0.94 


6.0 


2.2 


54.7 


b r 


s 


P 


n 


0.05 


22.4 


615. 


3.4 


0.94 


6.1 


2.3 


54.9 


b r 


s 


P 


n 


0.5 


21.8 


587. 


3.3 


0.95 


6.1 


1.9 


51.3 
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Chapter 7 

Comparison of All Flavors 



The new energy-exchange channels studied in the previous section lower the 
average energies. In order to compare the fluxes and spectra with those 
of v e and z/ e we perform a new series of runs where we include the full set of 
relevant microphysics for and also simulate the transport of z/ e and v e . 

Our findings for z/ e and z/ e agree with the previous literature, whereas our 
results for show lower mean energies. The mean energies of v e and are 
almost equal. The luminosity can differ from that of u e and v e by up to a 
factor of 2. Weak magnetism causes a difference between the mean energies 
of and of a few per cent. 

We present an overview of the results obtained by other groups and find 
that the more recent simulations agree rather well on the differential effects in 
energies and luminosities of the different flavors. These findings as well as our 
results suggest that assumptions commonly made for calculating oscillation 
effects, i.e., equal luminosities and a difference in the mean energies of v e and 

of up to a factor of 2, need to be changed. 

We show results of the detailed spectral shapes obtained in our simula- 
tions. Our a fit describes the spectra slightly better than a nominal Fermi- 
Dirac distribution. 

Parts of this chapter were already published in a similar form as: 
M. T. Keil, G. G. Raffelt, and H. T. Janka, "Monte Carlo study of supernova 
neutrino spectra formation," Astrophys. J. in press, astro-ph/0208035, 

7.1 Results from Our Monte Carlo Study 

For our v e and 9 e transport we employ the same microphysics as in Janka & 
Hillebrandt (1989a,b), i.e., charged- current reactions of with nucleons, iso- 
energetic scattering on nucleons, scattering on e 1 * 1 , and e + e~ pair annihilation. 
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In principle one should also include nucleon bremsstrahlung and the effect of 
nucleon recoils for the transport of u e and u e , but their effects will be minimal. 
Therefore, we preferred to leave the original working code unmodified for 
these flavors. 

In the first three rows of Table 17.11 we give the spectral characteristics 
for the Accretion-Phase Model I from the original simulation of Messer. The 
usual hierarchy of average neutrino energies is found, i.e., (e„ e ) : (e Pc ) : (e v ) = 
0.86 : 1 : 1.20. The luminosities are essentially equal between z/ e and z/ e while 
Vfx, v^i v T , and v T each provide about half of the v e luminosity. 

Our Monte Carlo runs of this profile establish the same picture for the 
same input physics (label "sp"). Although our mean energies are slightly 
offset to lower values for all flavors relative to the original run, our energies 
relative to each other are (e„ p ) : (e Pc ) : (e„ ) = 0.84 : 1 : 1.19 and thus very 
similar. However, once we include all energy exchanging processes (brspn) 
we find 0.84 : 1 : 1.02 instead. Therefore, (e v ) no longer exceeds (e Pe ) by 
much. The luminosity of is about half that of v c or z/ e which are approxi- 
mately equal, in rough agreement with the original results. Even though the 
additional processes lower the mean energy of they yield a more than 10% 
higher luminosity, mainly due to v e v e annihilation. 

Different rates for neutrinos and anti-neutrinos arise due to the weak 
magnetism correction of the neutrino-nucleon scattering cross section (Sec- 
tion 13. 3|) . The next two lines (wm) show results from runs where we included 
weak magnetism in our transport. Weak magnetism causes a significant 
correction to the neutrino-nucleon cross section that mainly arises due to 
the large anomalous magnetic moments of protons and neutrons (Vogel & 
Beacom 1999, Horowitz 2002). It increases the neutrino interaction rate but 
lowers the rate for anti-neutrinos. It is expected to be a small correction in 
the SN context, but has never been implemented for and so far. Fol- 
lowing Horowitz (2002) we add weak magnetism to our nucleon-recoil rate 
as given in Section EP1 

Our Monte Carlo code transports only one species of neutrinos at a time. 
In order to test the impact of weak magnetism we assumed that a chemical 
potential for would build up, and assumed a fixed value for the degen- 
eracy parameter throughout our stellar model. We then iterated several runs 
for and z/ M with different degeneracy parameters until their particle fluxes 
were equal because in a stationary state there will be no net flux of /z-lepton 
number. 

We performed this procedure for our Accretion-Phase Model I. The mean 
energy of v ^ goes down by 1% and goes up for by 4%. The mean lumi- 
nosities change in opposite directions by about 2-3%. Since we adjusted the 
particle fluxes by hand in order to fix the chemical potential, we cannot claim 
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Table 7.1. Comparing Monte Carlo results for different flavors. 
Accretion-Phase Models 



Flavor (e) flux (e 2 ) flux - — a p flux T n L 

\£p c /flux 



Accretion-Phase Model I 



Original 




















17.5 


388. 


1.20 


2.7 


0.97 


5.2 


1.1 


14.4 


v c 


14.6 


253. 


1 


4.4 


0.91 


3.5 


3.4 


29.2 




12.5 


190. 


0.86 


3.6 


0.93 


3.2 


2.8 


30.8 


Our runs 


















"n, v n ( S P) 


16.6 


362. 


1.19 


2.2 


1.01 


5.3 


-0.3 


15.8 


(brspn) 


14.3 


260. 


1.02 


2.7 


0.97 


4.3 


1.2 


17.9 


z/ M (wm) 


14.2 


254. 


1.01 


2.9 


0.96 


4.1 


1.6 


17.4 


z/ M (wm) 


14.9 


281. 


1.06 


2.8 


0.97 


4.3 


1.5 


18.3 




14.0 


237. 


1 


3.8 


0.93 


3.6 


2.7 


31.7 




11.8 


175. 


0.84 


2.9 


0.97 


3.4 


1.4 


31.9 


Accretion-Phase Model II 


Original 




















17.2 


380. 


1.09 


2.5 


0.98 


5.2 


0.8 


32.4 




15.8 


300. 


1 


4.0 


0.92 


4.0 


3.0 


68.1 




12.9 


207. 


0.82 


3.1 


0.96 


3.7 


1.7 


65.6 


Our runs 




















15.7 


317. 


1.02 


2.5 


0.98 


4.8 


0.8 


27.8 


Ve 


15.4 


283. 


1 


4.2 


0.92 


3.8 


3.2 


73.5 




13.0 


207. 


0.84 


3.4 


0.95 


3.6 


2.1 


73.9 


Note. — We give 


(e)flux 


and T in 


MeV, (e 2 ) 


flux in 


MeV 2 


, L v in 


10 51 erg s l . 
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20 40 60 80 100 
Radius [km] 



Fig. 7.1. — Comparison of the Accretion-Phase Model I calculations. Con- 
tinuous lines show our Monte Carlo runs while crosses represent the original 
simulation by Messer. The steps correspond to (e) = 3.15 T. Thin lines in 
the lower panel represent our findings after including the full set of relevant 
interaction processes. 

an effect on the luminosities. 

In Figure lTD we compare our calculations for the Accretion-Phase Model I 
with those of the original simulation. The step-like curve represents the mean 
energy of neutrinos in LTE for zero chemical potential. The smooth solid line 
is the mean energy (e) from our runs, the dotted line gives (e)fl ux . The crosses 
are the corresponding results from the original runs. In the case of (lower 
panel) we show the results of our "sp" run as thick lines and the results of our 
"brspn" run, i.e., the complete set of interactions included, as thin lines. The 
boundaries of our simulation domain were the same for all flavors, namely 
i?i n = 20 km and R OVLt = 100 km. 

As another example of an accreting proto-neutron star we use the 
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Accretion-Phase Model II. The neutrino interactions included in this model 
were nucleon bremsstrahlung, scattering on e 1 * 1 , and e + e~ annihilation. Nu- 
cleon correlations, effective mass, and recoil were taken into account, fol- 
lowing Burrows & Sawyer (1998, 1999), as well as weak magnetism effects 
(Horowitz 2002) and quenching of Qa at high densities (Carter & Prakash 
2002). All these improvements to the traditional microphysics affect mainly 

and to some degree also v e . Weak magnetism terms decrease the nucleon 
scattering cross sections for more strongly than they modify scatterings. 
In this hydrodynamic calculation, however, z/ M and were treated identically 
by using the average of the corresponding reaction cross sections. The effects 
of weak magnetism on the transport of and are therefore not included 
to very high accuracy. Note, moreover, that the original data come from a 
general relativistic hydrodynamic simulation with the solution of the Boltz- 
mann equation for neutrino transport calculated in the comoving frame of 
the stellar fluid. Therefore the neutrino results are affected by gravitational 
redshift and, depending on where they are measured, may also be blueshifted 
by Doppler effects due to the accretion flow to the nascent neutron star. 

Our Monte Carlo simulation in contrast was performed on a static back- 
ground without general relativistic corrections. It includes bremsstrahlung, 
recoil, e + e~ pair annihilation, scattering on e ± , and v e v e annihilation, i.e., our 
microphysics is similar but not identical with that used in the original run. 
As an outer radius we took 100 km; all flux parameters are measured at this 
radius because farther out Doppler effects of the original model would make 
it difficult to compare the results. Keeping in mind that we use very different 
numerical approaches and somewhat different input physics, the agreement 
in particular for v e and v e is remarkably good. This agreement shows once 
more that our Monte Carlo approach likely captures at least the differential 
effects of the new microphysics in a satisfactory manner. 

For the Accretion-Phase Model II we show the comparison with our runs 
in Figure 17.21 In this case, for the transport of our inner boundary is 
-R in = 16 km, while for v e and v e we use R in = 24 km. For v e and v e the 
charged-current processes (urea) keep these neutrinos in LTE up to larger 
radii than pair processes in the case of v^. With our choice of R in the neutrinos 
are in LTE within the innermost radial zones. In this profile our choice of 
boundaries reduces the CPU time needed and does not affect the results. As 
stated before outer boundaries were chosen to be at R out = 100 km. 

The results are similar to the Accretion-Phase Model I. The luminosities 
are not equipartitioned but instead follow roughly L Uc « L Pc « 2L^ M . The 
ratios of mean energies are (e„ e ) : (e Pe ) : (e„ ) = 0.82 : 1 : 1.09 in the original 
run and 0.84 : 1 : 1.02 in our run. 

In summary, both accretion-phase models agree reasonably well in the 
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Radius [km] 

Fig. 7.2. — Same as Figure I7~T1 for the Accretion-Phase Model II. 

(tve) '■ ( e P c ) ratio for all runs. Moreover, using traditional input physics one 
finds something like (ep c ) : (e u ) = 1 : 1.20. Depending on the implementa- 
tion of the new input physics and depending on the model one finds results 
between (e Pc ) : (e^) = 1 : 1.02 and 1 : 1.09. 

In order to estimate the corresponding results for later stages of the proto- 
neutron star evolution we employ our steep power-law model. We vary the 
power q of the temperature profile within a reasonable range so that q/p = 
0.25-0.35, with q and p defined in Equation (|4.5|) . Y e is fixed by demanding 
roughly equal number fluxes for u e and z/ e because a few seconds after bounce 
deleptonization should be essentially complete. The fluxes of these neutrinos 
depend very sensitively on Y e so that this constraint is only reached to within 
about 30% without tuning Y e to three decimal places. However, the mean 
energies are rather insensitive to the exact value of Y e . This is illustrated 
in the first section of Table I7~2*l bv the steep power-law model with q = 2.5 
where we show results for Y c = 0.15 and 0.20. The number fluxes of v e and v e 
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Table 7.2. Comparing Monte Carlo results for different flavors. 

Power-Law Models 



Flavor Y c (e) flux (e 2 ) flux - — ^- a p flux T r] L u 

\ £ pjflux 

Steep Power Law p = 10 



q = 2.5 







0.15 


20.4 


525. 


1.10 


2.8 


0.96 


5.9 


1.5 


23.5 




Ve 


0.15 


18.5 


413. 


f 


3.8 


0.92 


4.6 


3.0 


23.5 






0.15 


12.7 


198. 


0.69 


3.4 


0.94 


3.4 


2.4 


12.8 






0.2 


20.4 


521. 


1.14 


3.0 


0.97 


5.9 


1.5 


23.3 






0.2 


17.9 


383. 


f 


4.1 


0.92 


4.4 


3.1 


11.7 




V e 


0.2 


13.4 


218. 


0.75 


3.7 


0.93 


3.4 


2.9 


24.4 


q 


= 3.0 
























O.f 


17.7 


393. 


1.14 


2.9 


0.96 


5.0 


1.8 


12.7 




Ve 


O.f 


15.5 


289. 


f 


3.9 


0.93 


4.0 


2.8 


8.8 






O.f 


10.5 


132. 


0.68 


4.1 


0.92 


2.6 


3.0 


6.6 


q 


= 3.5 






















^, ^ 


0.07 


15.8 


310. 


1.22 


3.1 


0.95 


4.4 


2.1 


7.9 




z/ M (wm) 


0.07 


15.5 


296. 


1.19 


3.3 


0.94 


4.2 


2.3 


7.7 




z/ M (wm) 


0.07 


16.5 


337. 


1.27 


3.2 


0.95 


4.5 


2.1 


8.3 






0.07 


13.0 


207. 


f 


3.4 


0.94 


3.5 


2.3 


4.3 






0.07 


9.4 


103. 


0.72 


5.0 


0.90 


2.1 


3.9 


4.1 






Shallow Power Law 


p = 5, 


q = 


1 












0.3 


22.0 


596. 


1.14 


3.3 


0.94 


6.0 


2.2 


53.9 






0.3 


19.3 


440. 


f 


4.5 


0.91 


4.5 


3.7 


85.7 






0.3 


14.7 


262. 


0.76 


3.7 


0.93 


3.8 


2.7 


56.5 
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differ by less than 30% for Y c = 0.15, but differ by a factor of 3 for Y c = 0.2. 
At the same time, the average spectral energies barely change. 

The ratios of mean energies are not very different from those of the 
accretion-phase models. Of course, the absolute flux energies have no physi- 
cal meaning because we adjusted the stellar profile in order to obtain realistic 
values. For the luminosities we find L Uc < L v , different from the accretion 
phase. The steep power law implies that the radiating surfaces are similar for 
all flavors so that it is not surprising that the flavor with the largest energies 
also produces the largest luminosity. 

For our steepest profile we also performed runs including weak magnetism, 
because this should probe how much spectra are affected in late-time profiles. 
As already found for the Accretion-Phase Model I, the mean energy of goes 
up by 4%, whereas in the case of it decreases by 1%. The mean luminosities 
are almost unaffected. We conclude that weak-magnetism corrections are 
small. Transporting v^, and separately in a self-consistent hydrodynamic 
simulation is probably not worth the cost in computer time. 

We find that (e„ ) always exceeds (e Pe ) by a small amount, the exact value 
depending on the stellar model. During the accretion phase the energies seem 
to be almost identical, later they may differ by up to 20%. We have not found 
a model where the energies differ by the large amounts which are sometimes 
assumed in the literature. At late times when Y c is small the microphysics 
governing P e transport is closer to that for than at early times. Therefore, 
one expects that at late times the behavior of u c is more similar to than at 
early times. We do not see any argument for expecting an extreme hierarchy 
of energies at late times for self-consistent stellar models. 

We never find exact equipartition of the flavor- dependent luminosities. 
Depending on the stellar profile the fluxes can mutually differ by up to a 
factor of 2 in either direction. 

7.2 Previous Literature 

Studies of oscillation effects in SN neutrino spectra usually assume exactly 
equal neutrino luminosities for all flavors and a strong hierarchy for the mean 
energies. Often times they assume the difference between the mean energies 
of P e and to be about a factor of 2. We find an almost orthogonal picture 
in our simulations. Where does this come from? 

To the best of our knowledge, except in the very recent models by Buras et 
al. (2003b), the microphysics employed for transport was roughly the same 
in all published simulations. It included iso-energetic scattering on nucleons, 
e + e~ annihilation and scattering. Of course, the transport method and 



7.2. Previous Literature 



75 



the numerical implementation of the neutrino processes differ in the codes 
of different groups. The new results of self-consistent calculations point in 
the same direction as the findings from our Monte Carlo studies presented 
in the previous chapter. Therefore these findings are almost orthogonal to 
the previous assumptions oscillation studies are based on. Our representative 
sample of pertinent results is summarized in Table 1731 We have also included 
our Accretion-Phase Models I and II and a very recent model by the Garching 
group including all relevant processes. Note that the simulations discussed 
below did not in all cases use the same stellar models and equations of state 
for the dense matter in the SN core. 

We begin with the simulations of the Livermore group who find robust 
explosions by virtue of the neutron-finger convection phenomenon. Neutrino 
transport is treated in the hydrodynamic models with a multigroup flux- 
limited diffusion scheme. Mayle, Wilson, & Schramm (1987) gave detailed 
results for their SN simulation of a 25 M & star. For half a second after 
bounce they obtained a somewhat oscillatory behavior of the neutrino lu- 
minosities. After the prompt peak of the electron neutrino luminosity, they 
got L Ve w L Pc ps 2L;, M ps 50-130 x 10 51 erg s _1 . After about one second the 
values stabilize. This calculation did not produce the "standard" hierarchy 
of energies. However, there is clearly a tendency that 9 e behave more similar 
to at late times. 

The most recent published Livermore simulation is a 20 M star (Totani 
et al. 1998). It shows an astonishing degree of luminosity equipartition from 
the accretion phase throughout the early Kelvin-Helmholtz cooling phase. 
About two seconds after bounce the flux falls off more slowly than the 
other flavors. In Table 17.31 we show representative results for an early and a 
late time. The mean energies and their ratios are consistent with what we 
would have expected on the basis of our study. 

With a different numerical code, Bruenn (1987) found for a 25 M Q pro- 
genitor qualitatively different results for luminosities and energies. At about 
0.5 s after bounce the luminosities and energies became stable at the values 
given in Table 17.31 This simulation is an example for an extreme hierarchy 
of mean energies. 

In Burrows (1988) all luminosities are said to be equal. In addition it is 
stated that for the first 5 seconds (e u ) ~ 24 MeV and the relation to the 
other flavors is (e„ e ) : (e Pe ) : (e„ ) = 0.9 : 1 : 1.8. Detailed results are only 
given for u e , so we are not able to add this reference to our table. The large 
variety of models investigated by Burrows (1988) and the detailed results 
for v c go beyond the scope of our brief description. In a later paper Myra 
& Burrows (1990) studied a 13 M Q progenitor model and found the extreme 
hierarchy of energies shown in our table. 



Table 7.3. Flavor dependent flux characteristics from the literature. 
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Totani et al. (1998) 
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Bruenn (1987) 
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16 


Myra & Burrows (1990) 


0.13 
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Janka & Hillebrandt (1989b) 
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Accretion-Phase Model I (original) 0.32 

Accretion-Phase Model I (our run) 0.32 

Accretion-Phase Model II (original) 0.15 

Accretion-Phase Model II (our run) 0.15 

Buras et al. (personal comm.) 0.25 
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Table 7.3 — Continued 
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Mezzacappa et al. (2001) 0.5 


16 19 24 0.84 : 1 : 1.26 


25 


25 8 


Liebendorfer et al. (2001) 0.5 


19 21 24 0.90 : 1 : 1.14 


30 


30 10 



Note. — We give the time post bounce (tpb) in s, (e) in MeV, and L v in 10 51 erg s 
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With the original version of our code Janka & Hillebrandt (1989b) per- 
formed their analyses for a 20 M progenitor from a core-collapse calculation 
by Hillebrandt (1987). Of course, like our present study, these were Monte 
Carlo simulations on a fixed background model, not self-consistent simula- 
tions. Taking into account the different microphysics the mean energies are 
consistent with our present work. The mean energies of v e were somewhat on 
the low side relative to u e and the u e luminosity was overestimated. Both can 
be understood by the fact that the stellar background contained an overly 
large abundance of neutrons, because the model resulted from a post-bounce 
calculation which only included electron neutrino transport. 

Suzuki (1989) studied models with initial temperature and density pro- 
files typical of proto-neutron stars at the beginning of the Kelvin-Helmholtz 
cooling phase about half a second after bounce. He used the relatively stiff 
nuclear equation of state developed by Hillebrandt & Wolff (1985). In our 
table we show the results of the model C12. From Suzuki (1991) we took 
the model labeled C20 which includes bremsstrahlung. The model C48 from 
Suzuki (1993) includes multiple-scattering suppression of bremsstrahlung. 
Suzuki's models are the only ones from the previous literature which go be- 
yond the traditional microphysics for transport. It is reassuring that his 
ratios of mean energies come closest to the ones we find. 

Over the past few years, first results from Boltzmann solvers coupled 
with hydrodynamic simulations have become available, for example the un- 
published ones that we used as our Accretion-Phase Models I and II. For 
convenience we include them in Table 17.31 Further, we include a very recent 
accretion-phase model of the Garching group (Buras et al., personal commu- 
nication) that includes the full set of microphysical input. Finally, we include 
two simulations similar to the Accretion-Phase Model I, one by Mezzacappa 
et al. (2001) and the other by Liebendorfer et al. (2001). These latter papers 
show rms energies instead of mean energies. Recalling that the former tend to 
be about 45% larger than the latter these results are entirely consistent with 
our accretion-phase models. Moreover, the ratios of (e) rms tend to exaggerate 
the spread between the flavor- dependent mean energies because of different 
amounts of spectral pinching, i.e., different effective degeneracy parameters. 
To illustrate this point we take the first two rows from Table 17.11 as an ex- 
ample. The ratio of the mean energies (ei) = 17.5 MeV and (62) = 14.6 MeV 
is (ei)/(e 2 ) = 1.19. Using Equation 15. 141 for a\ = 2.7 and ct 2 = 4.4 the ratio 
of rms energies equals 1.31. 

To summarize, the frequently assumed exact equipartition of the emitted 
energy among all flavors appears only in some simulations of the Livermore 
group and some by Suzuki. We note that the flavor-dependent luminosities 
tend to be quite sensitive to the detailed atmospheric structure and chemical 
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composition. On the other hand, the often-assumed extreme hierarchy of 
mean energies was only found in the early simulations of Bruenn (1987) 
and of Myra & Burrows (1990), possibly a consequence of the neutron-star 
equation of state used in these calculations. 

If we ignore results which appear to be "outliers" , the picture emerging 
from Table 17.31 is quite consistent with our own findings. For the luminosi- 
ties, typically L Uc « L Dc and a factor of 2-3 between this and L Vfl in either 
direction, depending on the evolutionary phase. For the mean energies we 
read typical ratios in the range of (e„ e ) : (ep e ) : (e v ) = 0.8-0.9 : 1 : 1.0-1.3. 
The more recent simulations involving a Boltzmann solvers show a consis- 
tent behavior and will in future provide reliable information about neutrino 
fluxes and spectra. There is a clear tendency that the inclusion of all relevant 
microphysics decreases the hierarchy of mean energies. 

Analyzing detected neutrino spectra from a future galactic SN will de- 
mand very accurate predictions from self-consistent simulations. This is nec- 
essary for both, a possible measurement of neutrino oscillation parameters 
and also for testing the current SN paradigm. 

7.3 Spectral Shape 

Thus far we have characterized the neutrino spectra by a few simple parame- 
ters. However, it is extremely useful to have a simple analytic fit to the overall 
spectrum that can be used, for example, to simulate the response of a neu- 
trino detector to a SN signal. To study the quality of different fit functions 
we have performed a few high-statistics Monte Carlo runs for the Accretion- 
Phase Model I and two of the steep power-law models (p = 10, q — 2.5 and 
3.5) including all interaction processes. Moreover, we have performed these 
runs for the flavors v e , v e , and v^. 

In order to get smooth spectral curves we have averaged the output of 
70,000 time steps. In addition we have refined the energy grid of the neutrino 
interaction rates. Both measures leave the previous results unaffected but 
increase computing time and demand for memory significantly. 

In Figure 17.31 we show our high-statistics Monte Carlo spectra (MC) for 
the Accretion-Phase Model I together with the a-fit function f a (e) defined in 
Equation ()5.9|) and the //-fit function f v (e) of Equation (|5.15|) . The analytic 
functions can only fit the spectrum well over a certain range of energies. We 
have chosen to optimize the fit for the event spectrum in a detector, assuming 
the cross section scales with e 2 . Therefore, we actually show the neutrino flux 
spectra multiplied with e 2 . Accordingly, the parameters a and e, as well as r\ 
and T and the normalizations are determined such that the energy moments 
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(e 2 ), (e 3 ), and (e 4 ) are reproduced by the fits. 

Below each spectrum we show the ratio of our Monte Carlo results with 
the fit functions. In the energy range where the statistics in a detector would 
be reasonable for a galactic SN, say from 5-10 MeV up to around 40 MeV, 
both types of fits represent the Monte Carlo results nicely. However, in all 
cases the en-fit works somewhat better than the 77-fit. 

We have repeated this exercise for the steep power-law models with q = 
2.5 and the one with q = 3.5 and show the results in Figures 17.51 and \7A\ 
respectively. The quality of the fits is comparable to the previous example. 

7.4 Summary 

In comparing the spectra of all neutrino flavors we find a consistent behavior 
of all the considered models. Once all relevant interactions are included the 
differences between u e and spectra become rather small. Concerning mean 
energies from our two realistic profiles that represent the accretion phase, 
we always find (e„ e ) : (e Pe ) : (e„ M ) = 0.8-0.9 : 1 : 1.0-1.09. These findings 
are completely consistent with full hydrodynamic simulations that include 
the full set of interactions, i.e., the recent Garching-group models. For the 
luminosities we find good agreement, too. The runs commonly show L Uc w 
L Pc m 2 L Vfl . 

For our power-law models we find somewhat stronger hierarchies in the 
mean energies. It is not clear if this will still hold in self-consistent simula- 
tions. There are currently no self-consistent models with an accurate treat- 
ment of the neutrino transport available that continued the calculations to 
the late phases where steep proto-neutron-star profiles are found. 

In our collection of previous literature we found a variety of different 
predictions. The strong hierarchies in mean energies that were commonly 
assumed in the literature for calculating neutrino-oscillation effects are rather 
outliers. Also the exact equipartition of luminosities is not found in elaborate 
neutrino transport models. 

The standard picture of flavor-dependent SN neutrino spectra needs to be 
changed. Oscillation studies should allow for large differences in the v e and 

luminosities, and for almost equal mean energies of these flavors. Results 
of different self-consistent simulations show large variations in their absolute 
values. Therefore, only concepts that rely on differential effects of flavor- 
dependent SN neutrino spectra can have predictive power for identifying 
oscillation signatures. 
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Fig. 7.3. — High-statistics spectra for Accretion-Phase Model I including all 
interaction processes. The Monte Carlo (MC) results are shown as crosses, 
the analytic fit functions as smooth lines. The left-hand panels use as fits 
f a (e) according to Equation ()5.9jh the right-hand panels f v (e) according to 
Equation (|5.15|) . Below the spectra we show the ratio between Monte Carlo 
and fit. 
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Fig. 7.4. — Same as Figure [77U for the p — 10, q — 3.5 power law. 
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Fig. 7.5. — Same as Figure 17731 for the p — 10, q — 2.5 power law. 
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Chapter 8 



Detecting Oscillations of SN 
Neutrinos 

By now, neutrino oscillations are firmly established and the next goal is to pin 
down the oscillation parameters that are still unknown, like the neutrino mass 
hierarchy and the small mixing angle #13. Taking our findings on neutrino 
spectra into account requires new methods for extracting information, but 
enables us to use ratios of observables instead of their absolute magnitude. 

We present a method for identifying the effect of earth-matter oscillations 
in a detected SN neutrino signal by comparing the time dependence of the 
energy deposited in two detectors. By identifying spectral modulations with 
a Fourier-transform method the detection at a single detector is possible, 
but demands better energy resolution than the two-detector method. Finding 
oscillation signatures in the signal can pin down the mass hierarchy or restrict 
the small mixing angle. 

Parts of this chapter were published in a similar form in our publications 
A. S. Dighe, M. T. Keil, and G. G. Raffelt, 

"Detecting the neutrino mass hierarchy with a supernova at IceCube," JCAP, 
in press (2003), |hep-ph /0303210, and "Identifying earth matter effects on 
supernova neutrinos at a single detector," hep-ph/0304150, 

8.1 Oscillations of SN Neutrinos 

More and more accurate measurements have firmly established neutrino os- 
cillations. After many years of atmospheric- and solar-neutrino experiments 
and also long-baseline experiments (KamLAND and K2K) we have a good 
understanding of how neutrinos oscillate (Bahcall, Gonzales-Garcia, & Pena- 
Garay 2003, Fogli et al. 2002, Gonzales-Garcia & Nir 2002, de Holanda & 
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Smirnov 2003). The flavors u e , z/ M , and v T are in weak interaction eigenstates 
and therefore relevant when we talk about interactions of particles. Each of 
these weak eigenstates is a non-trivial superposition of three mass eigenstates 
z/i, z/ 2 , and u 3 , 



( Ue \ 








= u 




\ V rJ 







where U is the leptonic mixing matrix that can be written in the canonical 
form 

(I \ / ci 3 e i5 s l3 \ / c 12 8 X2 0\ 
U = c 23 s 23 \\ 1 -sx2 c l2 . (8.2) 

\0 -s 23 c 23 ) \-e~ iS s 13 c 13 / V 1/ 

Here Cx 2 = cos 6x2 and sx 2 = sin 9 12 etc., and 6 is a phase that can lead to 
CP- violating effects, that are, however, irrelevant for SN neutrinos. The three 
matrices then just correspond to the rotation matrices in the 2-3, 1-3, and 
1-2 planes, respectively. 

The mass eigenstates are most convenient for describing neutrino prop- 
agation in space. If there are three different mass eigenstates there can be 
three different masses. The oscillation, i.e., the interference of the three eigen- 
states, depends on the difference of the squares of their masses. Usually, the 
mass-squared differences are given as Am| = mf — m?. Of course, only two 
of them are independent. Compared to the current limit on the absolute 
neutrino mass, experiments find very small mass-squared differences. Two of 
the three masses appear to be rather degenerate. The third mass eigenstate 
can then lie far above or far below this doublet. Whether the third mass 
eigenstate lies above or below the almost degenerate doublet is referred to 
as normal or inverted hierarchy, respectively. Which hierarchy is realized in 
nature is an open question. 

The mass squared differences relevant for the atmospheric and solar neu- 
trino oscillations show the strong hierarchy Am^ tm ^> AmL This hierarchy, 
combined with the observed smallness of the angle #13 at CHOOZ (Apollonio 
et al. 1999) implies that the atmospheric neutrino oscillations essentially de- 
couple from the solar ones and each of these is dominated by only one of the 
mixing angles. The atmospheric neutrino oscillations are controlled by 9 23 
that may well be maximal (45°). The solar case is dominated by 8 12 , that is 
large but not maximal. To a reasonably good accuracy Am^ tm w Am^ and 
Am>Q m Awji and therefore Am| x m Am| 2 . From a global 3-flavor analysis 
of all data one finds the 3a ranges for the mass differences and mixing angles 
summarized in Table IH7T1 
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Obviously, in the context of SNe neutrino oscillations can only be found if 
flavor-dependent differences in the SN neutrino fluxes and spectra are present. 
We denote the fluxes of P e and at earth that would be observable in 
the absence of oscillations by F g ° and F^, respectively. In the presence of 
oscillations a v e detector actually observes 

F-F(e) = p D (e)F°(e) + [l - p D (e)} F° , (8.3) 

where p D (e) is the v e survival probability after propagation through the SN 
mantle and perhaps part of the earth before reaching the detector. Water 
Cherenkov detectors as well as scintillation detectors can only detect v e with 
a good efficiency. Therefore we give the z/ e flux at the detector. 

When neutrinos travel in matter their oscillatory behavior can change 
significantly. Depending on the neutrino source the path towards a detector 
leads through matter surrounding the source, like in stars or in a SN and 
might also go through the earth. This matter has a significant population of 
electrons, but other leptons are absent. By charged-current interactions the 
medium distinguishes between electron flavor and other flavors, which implies 
a potential for electron neutrinos and anti-neutrinos. In other words, the 
medium can change its refractive index for electron (anti-)neutrinos. Since the 
other flavors are unaffected, the oscillatory behavior can change significantly, 
depending on the electron density. 

A significant modification of the survival probability due to the propaga- 
tion through the earth appears only for those combinations of neutrino mixing 
parameters shown in Table 18.21 The earth matter effect depends strongly on 
two parameters, the sign of Am 2 2 and the value of |0i3| (Dighe & Smirnov 



Table 8.1. Neutrino mixing parameters 



Observation 


Mixing an 


gle 


Am 2 [meV 2 ] 


Sun, KamLAND 


012 = 27°- 


42° 


Am 2 21 = 


55-190 


Atmosphere, K2K 


023 = 32°- 


60° 


\Am 2 32 \ -- 


= 1400-6000 


CHOOZ 


013 < 14° 






Am 2 2 



Note. — Neutrino mixing parameters from a global anal- 
ysis of all experiments (3<r ranges) by Gonzalez- Garcia & 
Nir (2002). 
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2000, Dighe 2001). The normal hierarchy corresponds to m x < m 2 < m 3 , i.e., 
A77132 > 0, whereas the inverted hierarchy corresponds to 7713 < m\ < m-i-, 
i.e., Am| 2 < 0. Note that the presence or absence of the earth effect dis- 
criminates between values of sin 2 #13 less or greater than 10~ 3 , i.e., 6*13 less 
or larger than about 1.8°. Thus, the earth effect is sensitive to values of #13 
that are much smaller than the current limit fTable l8~T|) . 

Let us consider those scenarios where the mass hierarchy and the value of 
#13 are such that the earth effect appears for v e . In such cases the u e survival 
probability p D (E) is given by 



fP « cos 2 #12 - sin 29% sin(20® - 29 



12 



. 2 / Ami L 12.5 MeV\ , AS 

x sm 2 — , (8.4) 

V10" 5 eV 2 1000 km e ) y ' 

where the energy dependence of all quantities will always be implicit. Here 
9% is the mixing angle between v e and z/ 2 in earth matter while Am^ is the 
mass squared difference between the two anti-neutrino mass eigenstates V\ 
and T>2i L is the distance traveled through the earth, and e is the neutrino 
energy. We have assumed a constant matter density inside the earth, which 
is a good approximation for L < 10, 500 km, i.e., as long as the neutrinos do 
not pass through the core of the earth. 



Table 8.2. The earth effect in a SN signal. 



13- Mixing 


Normal Hierarchy 


Inverted Hierarchy 


sin 2 9 13 < 10- 3 


u e and v e 


u e and u e 


sin 2 9 n > 1Q- 3 







Note. — The earth effect appears for the indicated fla- 
vors in a SN signal. 
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8.2 Detecting Oscillations With Two 
Distant Detectors 

8.2.1 The Basic Principle 

How significant the oscillation signal in a detected SN neutrino spectrum 
is, basically depends on the difference between the u e and fluxes and the 
distance traveled through the earth. There is no exact prediction for the spec- 
tra emitted by a SN, but some features are established. The mean energies 
of v e and will be very similar, a difference of typically 0-20% should be 
expected. At the same time the number fluxes will be rather different by 
up to a factor 2 in any direction. The ratio of fluxes will certainly change 
with time, because there are two distinct phases responsible for the neu- 
trino signal, namely the accretion phase and Kelvin-Helmholtz cooling of the 
proto- neutron star (Chapter EJ). 

With the time dependence of the flux difference it is possible to find an 
oscillation signature without a detailed knowledge of the original SN signal. 
Oscillations inside the earth depend on the flux difference and will therefore 
also change with time. Comparing the time dependence of a neutrino signal 
of two detectors yields a good chance for detecting the earth effect. The setup 
should be as follows: one detector sees the SN from above, i.e., without earth 
effect, and another detector that neutrinos reach after traveling through the 
earth. The chance for encountering this setup increases with the separation 
of the detectors. Detectors that are separated by a long distance and will be 
operational for many decades are the future IceCube detector in Antarctica 
(Ahrens et al. 2002b) and, for example, the Super-Kamiokande detector in 
Japan. 

8.2.2 The SN Signal at IceCube 

It has been recognized for a long time that neutrino telescopes detecting 
Cherenkov light in ice can detect a SN neutrino burst because the Cherenkov 
glow of the ice can be identified as time-correlated noise among all phototubes 
(Halzen, Jacobsen, & Zas 1994, 1996). This approach has been used by the 
neutrino telescope AMANDA to exclude the occurrence of a galactic SN over 
a recent observation period (Ahrens et al. 2002a). 

The SN neutrinos streaming through the antarctic ice interact according 
to z/ c p — > ne + and some other less important reactions. The positrons, in 
turn, emit Cherenkov light producing a homogeneous and isotropic glow of 
the ice. The optical modules (OMs) that are frozen into the ice are immersed 
in this diffuse bath of photons and pick up a number corresponding to their 
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angular acceptance and quantum efficiency. Estimating the event rate of a 
SN signal induced by the z/ e p — > ne + reaction yields (Dighe et al. 2003a) 

r fio -i L » c /I0kpc\ 2 

1 events - 62 S 1Q 52 erg g -l [~1T~ ) ^ ' ^ } 

where L Dc is the u e luminosity after flavor oscillations, and D the distance 
between SN and detector. The factors /a ux and /d e t parameterize the flux 
and detector characteristics, respectively. Both factors are of order 1. For the 
flux we have 

15MeV 8 (e 3 c > 

Jflux — 



(e Pc ) 15 (15 MeV) 3 
(3 + a)(2 + a) 8 (e Pc > 2 



;i + a) 2 15 (15 MeV) 2 ' 



5.6) 



where we used Equation ()5.13|) in the last line, i.e., we applied our a param- 
eterization. For a = 3 and (ep c ) = 15 MeV the factor is 1. 

The factor /det contains the efficiency for converting the neutrino signal 
into photons by the ice as well as the efficiency for detecting these photons. 
For our purpose we assume /d e t = 1- 

In the 4800 OMs of IceCube a SN in our galaxy would yield a total event 
number of about 1.5 x 10 6 photons if we assume the SN radiates 5 x 10 52 erg. 
This rate needs to be compared to the background counting rate of 300 Hz 
per OM and thus for a duration of 10 s we expect 1.44 x 10 7 photons in total. 
Assuming Poisson fluctuations, the uncertainty of this number is 3.8 x 10 3 , 
i.e., 0.1% of the SN signal. Therefore, one can determine the SN signal with 
a statistical sub-percent precision, ignoring for now problems of absolute 
detector calibration. 

In order to illustrate the statistical power of IceCube to observe a SN 
signal we use two different numerical SN simulations. The first was performed 
by the Livermore group (Totani et al. 1998) that involves traditional input 
physics for interactions and a flux-limited diffusion scheme for treating 
neutrino transport. The great advantage of this simulation is that it covers 
the full evolution from infall over the explosion to the Kelvin-Helmholtz 
cooling phase of the proto-neutron star. We show the Livermore v e and z/ M 
lightcurves in Figure l8~T1 (left panels). For all flavors these curves are also 
displayed in Figure l2~4l 

Our second simulation was performed with the Garching code (Rampp 
& Janka 2002). It includes all relevant neutrino interaction rates, including 
nucleon bremsstrahlung, neutrino pair processes, weak magnetism, nucleon 
recoils, and nuclear correlation effects. The neutrino transport part is based 



8.2. Detecting Oscillations With Two Distant Detectors 



91 




Fig. 8.1. — Supernova u e and v ^ light curves and average energies. Left: Liv- 
ermore simulation (Totani et al.1999). Right: Garching simulation (Raffelt et 
al. 2003). 



on a Boltzmann solver. The neutrino-radiation hydrodynamics program al- 
lows one to perform spherically symmetric as well as multi-dimensional sim- 
ulations. The progenitor model is a 15 M Q star with a 1.28 M Q iron core. 
The period from shock formation to 468 ms after bounce was evolved in two 
dimensions. The subsequent evolution of the model is simulated in spherical 
symmetry. At 150 ms the explosion sets in, although a small modification 
of the Boltzmann transport was necessary to allow this to happen (Janka et 
al. 2003). Recall that unmanipulated full-scale models with an accurate treat- 
ment of the microphysics currently do not obtain explosions (Section I2.4j) . 
This run will be continued beyond the current epoch of 750 ms post bounce; 
we here use the preliminary results currently available (Raffelt et al. 2003). 
We show the Garching v e and lightcurves in Figure I8~T1 fright panels). 

We take the Livermore simulation to represent traditional predictions for 
flavor-dependent SN neutrino fluxes and spectra that were used in many 
previous discussions of SN neutrino oscillations. The Garching simulation is 
taken to represent a situation when the interactions are more system- 
atically included so that the flavor- dependent spectra and fluxes are more 
similar than had been assumed previously. We think it is useful to juxtapose 
the IceCube response for both cases. 

Another difference is that in Livermore the accretion phase lasts longer. 
Since the explosion mechanism is not finally settled, it is not obvious which 
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Fig. 8.2. — Supernova signal in IceCube assuming a distance of 10 kpc, based 
on the Livermore simulation (left) and the Garching one (right), in both cases 
ignoring flavor oscillations. In the bottom panels we have used 50 ms bins 
and have added noise from a background rate of 300 Hz per OM. 



case is more realistic. Moreover, there could be differences between different 
SNe. The overall features are certainly comparable between the two simula- 
tions. 

In Figure 18.21 we show the expected counting rates in IceCube as given 
in Equation (|8.5J) for an assumed distance of 10 kpc and 4800 OMs for the 
Livermore (left) and Garching (right) simulations. We also show this signal 
in 50 ms bins where we have added noise from a background of 300 Hz per 
OM. The baseline is at the average background rate so that negative counts 
correspond to downward background fluctuations. 

One could easily identify the existence and duration of the accretion phase 
and thus test the standard delayed-explosion scenario. One could also mea- 
sure the overall duration of the cooling phase and thus exclude the presence of 
significant exotic energy losses. Therefore, many of the particle-physics limits 
based on the SN 1987A neutrinos (Raffelt 1999) could be supported with a 
statistically serious signal. If the SN core were to collapse to a black hole after 
some time, the sudden turn-off of the neutrino flux could be identified. In 
short, when a galactic SN occurs, IceCube is a powerful stand-alone neutrino 
detector, providing us with a plethora of information that is of fundamental 
astrophysical and particle-physics interest. 
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8.2.3 The Oscillation Signal at Ice Cube 

In order to calculate the extent of the earth effect for IceCube, we will as- 
sume that the relevant mixing parameters are AmJ 2 = 6 x 10~ 5 eV 2 and 
sin 2 (2# 12 ) = 0.9. We further assume that the source spectra are given by the 
functional form of our a fit Equation (|5.9jl . The values of the parameters a 
and (e) for both the v e and spectra are in general time dependent. 

In Figure 18.31 we show the variation of the expected IceCube signal with 
earth-crossing length L for the two sets of parameters detailed in Table 18.31 
The first could be representative of the accretion phase, the second of the 
cooling signal. We use the two-density approximation for the earth density 
profile, where the core has a density of 11.5 g cm -3 and a radius of 3,500 km, 
while the density of the earth mantle was taken to be 4.5 g cm -3 . We observe 
that for short distances, corresponding to near-horizontal neutrino trajecto- 
ries, the signal varies strongly with L. Between about 3,000 and 10,500 km 
it reaches an asymptotic value that we call the "asymptotic mantle value." 
For Case (a), this value corresponds to about 1.5% depletion of the signal, 
whereas for (b) it corresponds to about 6.5% depletion. 

Beyond an earth-crossing length of about 10,500 km, the neutrinos have 
to cross the earth core with another large jump in density. The core ef- 




Earth— Crossing Length [km] 

Fig. 8.3. — Variation of the expected IceCube signal with neutrino earth 
crossing length L for the assumed flux and mixing parameters of Table 18.31 
The signal is normalized to 1 when no earth effect is present, i.e. for L = 0. 
The dashed line is for the case representing the accretion phase, the solid line 
for the cooling phase. 
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Fig. 8.4. — Asymptotic IceCube signal modification by the earth effect. The 
fixed flux parameters are (a) (e Pc ) = 15 MeV, a Pc = 4.0, and a p = 3.0 
and (b) (e Ce ) = 15 MeV, a Pc = a Ufi = 3.0. The contours are equally spaced 
starting from 1.02 (light) in 0.02 decrements to smaller values (darker). 

fects change the asymptotic mantle value by roughly 1% as can be seen in 
Figure 18.31 We neglect the core effects in the following analysis, and the 
"asymptotic value" always refers to the asymptotic mantle value. 

For the largest part of the sky the earth effect either appears with this 
asymptotic value ("neutrinos coming from below"), or it does not appear 
at all ("neutrinos from above"). Therefore, we now focus on the asymptotic 
value and study how the signal modification depends on the assumed flux 
parameters. In Table l8~4l we show the signal modification for (e Pe ) = 15 MeV, 
a Pc = 4.0, and a p ^ = 3.0 as a function of (e^) and the flux ratio $p c /$p^. In 
Table 15*31 we show the same with a Pc = a Pfi = 3.0. The results are shown in 
the form of contour plots in Figure 18.41 

Even for mildly different fluxes or spectra the signal modification is several 



Table 8.3. Flux parameters for two representative cases. 



Example 


Phase 




(%) 






$o / $ o 


Asymptotic 






[MeV] 


[MeV] 








Earth Effect 


(a) 


Accretion 


15 


17 


4 


3 


1.5 


-1.5% 


(b) 


Cooling 


15 


18 


3 


3 


0.8 


-6.5% 



8.2. Detecting Oscillations With Two Distant Detectors 



95 



Table 8.4. Asymptotic IceCube signal modification by the earth effect 



Flux ratio 






(%> 


[MeV] 








15 


16 


17 


18 


19 


20 


2.0 


1.026 


1.014 


1.002 


0.988 


0.974 


0.960 


1.9 


1.023 


1.011 


0.999 


0.985 


0.971 


0.956 


1.8 


1.021 


1.009 


0.995 


0.982 


0.967 


0.952 


1.7 


1.018 


1.005 


0.992 


0.978 


0.963 


0.948 


1.6 


1.015 


1.002 


0.988 


0.974 


0.959 


0.944 


1.5 


1.012 


0.998 


0.984 


0.969 


0.954 


0.939 


1.4 


1.008 


0.994 


0.980 


0.965 


0.949 


0.934 


1.3 


1.004 


0.990 


0.975 


0.960 


0.944 


0.928 


1.2 


1.000 


0.985 


0.970 


0.954 


0.938 


0.922 


1.1 


0.995 


0.980 


0.964 


0.948 


0.932 


0.915 


1.0 


0.989 


0.974 


0.957 


0.941 


0.925 


0.908 


0.9 


0.983 


0.967 


0.950 


0.934 


0.917 


0.901 


0.8 


0.976 


0.959 


0.942 


0.925 


0.909 


0.892 


0.7 


0.967 


0.950 


0.933 


0.916 


0.899 


0.883 


0.6 


0.958 


0.940 


0.923 


0.906 


0.889 


0.873 


0.5 


0.946 


0.928 


0.911 


0.894 


0.877 


0.862 



Note. — The fixed flux parameters are (e Pc ) = 15 MeV, 
a Pc = 4.0, and a Pfi = 3.0. 
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Table 8.5. Same as Table 18731 with a Pc = a p „ = 3.0. 



Flux ratio 




(%) [MeV] 




®V®1 15 


16 


17 18 19 


20 



2.0 


1 


.036 


1 


.024 


1 


.012 


1 


.000 





.986 





.972 


1.9 


1 


.033 


1 


.022 


1 


.010 





.996 





.983 





.968 


1.8 


1 


.031 


1 


.019 


1 


.006 





.993 





.979 





.964 


1.7 


1 


.028 


1 


.016 


1 


.003 





.989 





.975 





.960 


1.6 


1 


.025 


1 


.013 





.999 





.985 





.971 





.955 


1.5 


1 


.022 


1 


.009 





.995 





.981 





.966 





.951 


1.4 


1 


.019 


1 


.005 





.991 





.976 





.961 





.945 


1.3 


1 


.015 


1 


.001 





.986 





.971 





.955 





.940 


1.2 


1 


.010 





.996 





.981 





.965 





.949 





.933 


1.1 


1 


.006 





.991 





.975 





.959 





.943 





.927 


1.0 


1 


.000 





.985 





.969 





.952 





.936 





.919 


0.9 





.994 





.978 





.961 





.945 





.928 





.911 


0.8 





.986 





.970 





.953 





.936 





.919 





.903 


0.7 





.978 





.961 





.944 





.926 





.910 





.893 


0.6 





.968 





.950 





.933 





.916 





.899 





.882 


0.5 





.956 





.938 





.920 





.903 





.886 





.870 
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percent, by far exceeding the statistical uncertainty of the IceCube signal, 
although the absolute calibration of IceCube may remain uncertain to within 
several percent. However, the signal modification will vary with time during 
the SN burst. During the early accretion phase that is expected to last for 
a few 100 ms and corresponds to a significant fraction of the overall signal, 
the flux may be almost a factor of 2 smaller than the v e flux, but it will 
be slightly hotter and less pinched (Raffelt et al. 2003). This corresponds 
to Case (a) of Table 18. 3j it is evident from Figure 18.31 and Table 18.41 that 
this implies that the earth effect is very small. During the Kelvin-Helmholtz 
cooling phase the flux ratio is reversed with more being emitted than u e , 
but still with the same hierarchy of energies. This corresponds to Case (b); in 
this case the earth effect could be about 6%. This time dependence may allow 
one to detect the earth effect without a precise absolute detector calibration. 

In order to illustrate the time dependence of the earth effect we show 
in Figure IH31 the expected counting rate in IceCube for both the Livermore 
(left panels) and Garching (right panels) simulations. In the upper panels we 
show the expected counting rate with flavor oscillations in the SN mantle, 
but no earth effect (solid lines), or with the asymptotic earth effect (dashed 
lines) that obtains for a large earth-crossing path. Naturally the differences 
are very small so that we show in the lower panels the ratio of these curves, 
i.e., the expected counting rate with/without earth effect as a function of 
time for both Livermore and Garching. While for the Livermore simulation 
there is a large earth effect even at early times, the change from early to late 
times in both cases is around 4-5%. Therefore, the most model-independent 
signature is a time variation of the earth effect during the SN neutrino signal. 

In order to demonstrate the statistical significance of these effects we in- 
tegrate the expected signal for both simulations separately for the accretion 
phase and the subsequent cooling phase; the results are shown in Table ESI 
For both simulations the earth effect itself and its change with time is statis- 
tically highly significant. Due to more schematic input physics the Livermore 
model overestimates the differences between v e and spectra. However, the 
relative change of the earth effect during accretion and cooling is not vastly 
different between the two simulations. Recalling that the absolute detector 
calibration may be very uncertain so that one has to rely on the temporal 
variation of the earth effect, the difference between Livermore and Garching 
becomes much smaller. We expect that it is quite generic that the temporal 
change of the earth effect is a few percent of the overall counting rate. 
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Fig. 8.5. — Earth effect in IceCube. The upper panels show the expected 
counting rate based on the Livermore (left) and Garching (right) models, 
including flavor oscillations. The solid line is without earth effect [L = 0), 
the dashed line with asymptotic earth effect (L ^ 0). The lower panels show 
the ratio between these curves, i.e., the ratio of counting rates with/without 
earth effect. 



Table 8.6. IceCube Cherenkov counts for the numerical SN models. 



Livermore Garching 





Accretion 


Cooling 


Accretion 


Cooling 


Integration time [s] 


0-0.500 


0.500-3 


0-0.250 


0.250-0.700 


SN Signal [Counts] 










No Earth Effect 


519,080 


818,043 


173,085 


407,715 


Asymptotic Earth Effect 


488,093 


751,137 


171,310 


390,252 


Difference 


30,987 


66,906 


1,775 


17,463 


Fractional Difference 


-5.97% 


-8.18% 


-1.03% 


-4.28% 


Background [Counts] 


720,000 


4,320,000 


360,000 


648,000 


^/Background / Signal 


0.16% 


0.25% 


0.35% 


0.20% 
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8.2.4 Super- or Hyper-Kamiokande and IceCube 

One can measure the earth effect in IceCube only in conjunction with another 
high-statistics detector. We do not attempt to simulate in detail the SN sig- 
nal in this other detector but simply assume that it can be measured with a 
precision about as good as in IceCube. One candidate is Super-Kamiokande, 
a water Cherenkov detector that would measure around 10 4 events from a 
galactic SN at a distance of 10 kpc. Therefore, the statistical precision for the 
total neutrino energy deposition in the water is around 1% and thus worse 
than in IceCube. Even though Super-Kamiokande will measure a larger num- 
ber of Cherenkov photons than IceCube, a single neutrino event will cause 
an entire Cherenkov ring to be measured, i.e., the photons are highly cor- 
related. Therefore, in the estimated statistical \/N fluctuation of the signal, 
the fluctuating number N is that of the detected neutrinos. If the future 
Hyper-Kamiokande is built, its fiducial volume would be about 30 times that 
of Super-Kamiokande. In this case the statistical signal precision exceeds that 
of IceCube for the equivalent observable. 

We denote the equivalent IceCube signal measured by Super- or Hyper- 
Kamiokande as iVsK and the IceCube signal as N\c- If the distances traveled 
by the neutrinos before reaching these two detectors are different, the earth 
effect on the neutrino spectra may be different, which will reflect in the ratio 
N SK /Ni C . Of course, in the absence of the earth effect this ratio equals unity 
by definition. 

The geographical position of IceCube with respect to Super- or Hyper- 
Kamiokande at a northern latitude of 36.4° is well-suited for the detection 
of the earth effect through a combination of the signals. Using Figure 18.31 
we can already draw some qualitative conclusions about the ratio Nsk/Nic- 
Clearly, Nqk/Nic = 1 if neutrinos do not travel through the earth before 
reaching either detector. If the distance traveled by neutrinos through the 
earth is more that 3,000 km for both detectors, the earth effects on both N$& 
and iVic are nearly equal and their ratio stays around unity. If the neutrinos 
come "from above" for SK and "from below" for IceCube, or vice versa, the 
earth matter effect will shift this ratio from unity. 

In Figure EH we show contours of Nsk/Njc for the SN position in terms 
of the location on earth where the SN is at the zenith. The map is an area 
preserving Hammer-Aitoff projection so that the sizes of different regions in 
the figure gives a realistic idea of the "good" and "bad" regions of the sky. In 
order to generate the contours we use the parameters of Case (b) in Table 
so that the asymptotic suppression of the signal is about 6.5%. The sky falls 
into four distinct regions depending on the direction of the neutrinos relative 
to either detector as described in Table IH771 When the neutrinos come from 
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above for both detectors (Region D) there is no earth effect. If they come 
from below in both (Region C), the earth effect is large in both. Depending 
on the exact distance traveled through the earth, the event ratio can be large, 
but generally fluctuates around 1. In the other regions where the neutrinos 
come from above for one detector and from below for the other (Regions A 
and B) the relative effect is large. 




Table 8.7. Regions in Figure l8~T)l for the earth effect in IceCube and 

Super-Kamiokande. 



Region 


Sky fraction 


Neutrinos 


come from 


Nsk/Njc 






IceCube 


Super-K 




A 


0.35 


below 


above 


1.070 


B 


0.35 


above 


below 


0.935 


C 


0.15 


below 


below 


Fluctuations around 1 


D 


0.15 


above 


above 


1 
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8.3 Detecting the Earth- Matter Effect at a 
Single Detector 

According to Equation (|8.4j) the modulations imprinted on the detected spec- 
tra are proportional to sin 2 (#/e). Therefore the prefactor # should appear 
as a dominant frequency in the Fourier transform of the inverse energy spec- 
trum (Dighe et al. 2003b), because in inverse-energy space the modulations 
of the SN neutrino spectrum are nearly equispaced. 

The equidistant peaks in the modulation of the inverse-energy spectrum 
are a necessary feature of the earth effects. Indeed, the net v e flux at the 
detector may be written using ()8.3j) and (|8.4jl in the form 

FP = sin 2 12 F° + cos 2 6 12 F° + AF% sin 2 (Am|L/e) , (8.7) 

where AF° = (F E ° — Fj[) depends only on the primary neutrino spectra, 
whereas = — sin 29% sin(29® 2 — 2# 12 ) depends only on the mixing pa- 
rameters and is independent of the primary spectra. The last term in ()8.7|) 
is the earth oscillation term that contains a frequency fc e = 2Am 2 B L in e _1 , 
with the coefficient AF°A& being a comparatively slowly varying function of 
e _1 . The first two terms in (J8.7)) are also slowly varying functions of e _1 , and 
hence contain frequencies in e _1 that are much smaller than 

The frequency is completely independent of the primary neutrino spec- 
tra, and indeed can be determined to a good accuracy from the knowledge of 
the solar oscillation parameters, the earth matter density, and the direction 
of the SN. If this frequency component is isolated from the inverse-energy 
spectrum of u e , the earth effects would be identified. 

Due to the energy dependence of the prefactors of the oscillation term 
AF°A S) and of Ami in Equation ()8.7|) the peak in the Fourier-transformed 
spectrum gets a certain width around k®. This peak can be clearly identified 
if the earth effect is present. 

In an experiment there are two effects obscuring the signal. One is the 
statistical fluctuations of the signal and the other is the smearing of the 
modulation signal by the energy resolution of the detector. Higher statistics 
as well as better energy resolution enhance the oscillation signature. A pre- 
scription for identifying the prominent peak on top of the background can 
be found in Dighe et al. (2003b). If the direction of the SN is known we can 
estimate k^ and by that simplify the peak identification. 

Once the peak is identified this is a clear signal of earth-matter effects. 
As stated in Section 18.2.41 an identification of the earth-matter effect will 
severely restrict the neutrino mixing parameter space, since the effects are 
present only with the combinations of the neutrino mass hierarchy and the 
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mixing angle #13 given in Table 18.21 In particular, if sin 2 9 13 is measured at a 
laboratory experiment to be greater than 10 -3 , then the earth effects on the 
i> e spectrum imply the normal mass hierarchy. However, if the earth effects 
are not detected, it does not rule out any neutrino mixing parameters, owing 
to the current uncertainties in the primary fluxes. 

With the position of the peak it is even possible to determine AtUq to an 
accuracy of a few percent and therefore much better than the current limits, 
that are Awig = (5.5-19) x 10~ 5 eV 2 . The Fourier-transform method will 
measure AtUq to a precision of 10%, comparable to what can be reached by 
KamLAND (de Gouvea & Pena-Garay, 2001). 

8.4 Summary 

With our new findings for the SN neutrino spectra methods for identifying 
oscillation signatures in an observed neutrino spectrum need to be changed. 
We presented two very promising concepts for detecting the earth-matter ef- 
fect. For the IceCube detector in Antarctica earth matter-effects are present 
in the signal of a future galactic SN on the level of a few percent. If the Ice- 
Cube signal can be compared with another high-statistics signal, notably in 
Super-Kamiokande or Hyper-Kamiokande, the earth effect becomes clearly 
visible as a difference between the detectors. As one is looking for a signal 
modification in the range of a few percent, the absolute detector calibration 
may not be good enough in one or both of the instruments. However, for typ- 
ical numerical SN simulations the effect is time dependent and most notably 
differs between the early accretion phase and the subsequent neutron star 
cooling phase. Therefore, one would have to search for a temporal variation 
of the relative detector signals of a few percent. The large number of optical 
modules in IceCube renders this task statistically possible. In fact depend- 
ing on the differences in flavor-dependent fluxes, the statistical accuracy of 
Super-Kamiokande may turn out to be the limiting factor. This limitation is 
not significant for Hyper-Kamiokande. 

The unique location of IceCube in Antarctica implies that for about 70% 
of the sky this detector sees the SN through the earth when Super- and Hyper- 
Kamiokande see it from above, or the other way round, i.e., the chances of 
a relative signal difference between the detectors are large. If both detectors 
were to see the SN from above or both through the earth, the comparison of 
the signals would would not reveal the earth effect. 

Once future detectors like Hyper-Kamiokande or large scintillation de- 
tectors become available a powerful way of identifying earth-matter effects is 
through the Fourier transform of the inverse-energy spectrum. In the inverse- 
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energy spectrum modulations due to earth matter will be approximately eq- 
uispaced and therefore appear as a single peak in the Fourier transformed 
spectrum. With this method one can unambiguously prove the existence of 
earth-matter oscillations and even measure AtUq to a accuracy of a few per- 
cent. 

Identifying the earth-matter effects in the neutrino signal of a future 
galactic SN would restrict the parameter space for neutrino oscillations 
severely. If, in addition, the magnitude of the mixing angle #13 can be estab- 
lished to be large in the sense of sin 2 9 13 ~ 10~ 3 by a long-baseline experiment 
(Barger et al. 2001, Cervera et al. 2000, Freund et al. 2001), it implies the 
normal mass hierarchy. 

On the other hand, if sin 2 6> 13 < 10~ 3 has been established, the earth 
effect is unavoidable whatever the hierarchy is. Not observing it for such a 
small #13 would imply that the primary SN neutrino fluxes and spectra are 
more similar than indicated by state-of-the-art numerical simulations. For 
sin 2 #13 > 10~ 3 not observing the earth- matter effects does not allow one to 
draw conclusions, because it can be due the original SN spectra and fluxes, 
or due to the neutrino mass hierarchy. 
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Chapter 9 

Discussion and Summary 



It has been known for some time that the "traditional" set of i^-matter in- 
teractions employed by hydrodynamic SN simulations was incomplete. In 
addition to the ingredients that were known to be missing, i.e., nucleon 
bremsstrahlung and nucleon recoil, we showed that u c u c annihilation into 
v^Vfj, is always more important than the traditional e + e~ annihilation pro- 
cess by a factor of 2-3. 

In a systematic approach we studied the formation of neutrino spectra 
and fluxes in a SN core. Using a Monte Carlo code for neutrino transport, 
we varied the microscopic input physics as well as the underlying static 
proto-neutron star atmosphere. We used two background models from self- 
consistent hydrodynamic simulations, and several power-law models with 
varying power-law indices for the density and temperature and different con- 
stant values for the electron fraction Y e . 

The transport opacity is dominated by neutral-current scattering 
on nucleons. In addition, there are number-changing processes (nucleon 
bremsstrahlung, leptonic pair annihilation) and energy- changing processes 
and v^Ve, v^v e scattering). Recoil in nucleon scattering allows for a 
small energy exchange in each collision. The spectra and fluxes are roughly 
accounted for if one includes one significant channel of pair production and 
one for energy exchange in addition to z/^N scattering. For example, the tra- 
ditional set of microphysics (iso-energetic z/ M N scattering, e + e~ annihilation, 
and scattering) yields crudely comparable spectra and fluxes to a cal- 
culation where pairs are produced by nucleon bremsstrahlung and energy 
is exchanged by nucleon recoil. The overall result is robust to within 30% 
against the detailed choice of microphysics. 

However, in view of neutrino oscillations, where flavor- dependent flux 
differences are important, state-of-the-art simulations should aim at a high 
precision for the fluxes and spectral energies. Therefore, one needs to include 
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bremsstrahlung, leptonic pair annihilation, neutrino-electron scattering, and 
energy transfer in neutrino-nucleon collisions. As expected, the traditional 
e + e~ annihilation process is always much less important than v e v c annihila- 
tion. None of the reactions studied here can be neglected except perhaps the 
traditional e + e~ annihilation process, and u^u c and v^v^ scattering. 

The existing treatments of the nuclear-physics aspects of the NN — ► NNz/z/ 
bremsstrahlung process are rather schematic. We find, however, that the 

fluxes and spectra do not depend sensitively on the exact strength of 
the bremsstrahlung rate. Therefore, while a more adequate treatment of 
bremsstrahlung remains desirable, the final results are unlikely to be much 
affected. 

The transport of and is usually treated identically. However, weak- 
magnetism effects render the z/^N and z/^N scattering cross sections somewhat 
different (Horowitz 2002), causing a small chemical potential to build up. 
We find that the differences between the average energies of and are 
only a few percent and can thus be neglected for most purposes. 

Including all processes works in the direction of making the fluxes and 
spectra of more similar to those of v c compared to a calculation with the 
traditional set of input physics. During the accretion phase the neutron-star 
atmosphere is relatively expanded, i.e., the density and temperature gradients 
are relatively shallow. Our investigation suggests that during this phase (e v ) 
is only slightly larger than (e Po ), perhaps by a few percent or 10% at most. 
This result agrees with the first hydrodynamic simulation including all of 
the relevant microphysics except v e v e annihilation (Accretion-Phase Model 
II) provided to us by M. Rampp. For the luminosities of the different neutrino 
species one finds L Pc pa L Uc pa 2L Ufi . The smallness of L Vfi is not surprising 
because the effective radiating surface is much smaller than for P e . 

During the Kelvin-Helmholtz cooling phase the neutron-star atmosphere 
will be more compact, the density and temperature gradients will be steeper. 
Therefore, the radiating surfaces for all species will become more similar. 
In this situation L v may well become larger than L Pc . However, the rela- 
tive luminosities depend sensitively on the electron concentration. Therefore, 
without a self-consistent hydrostatic late-time model it is difficult to claim 
this luminosity cross-over with confidence. 

The ratio of the spectral energies is most sensitive to the temperature 
gradient relative to the density gradient. In our power-law models we used 
p oc r~ p and T oc r~ q . Varying q/p between 0.25 and 0.35 we find that (e Pe ) : 
(e v ) varies between 1 : 1.10 and 1 : 1.22. Noting that the upper range for q/p 
seems unrealistically large we conclude that even at late times the spectral 
differences should be small; 20% sounds like a safe upper limit. However, the 
power-law models might overestimate the spread of mean energies, as can 
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be inferred from comparing the shallow power-law model with the realistic 
accretion-phase models. We stress that there is no physical reason for the 
inequality (e Pc ) < {e v ). It may well be that at some point the mean energies 
are equal or even cross over. 

The statements about the relative neutrino energies in the previous lit- 
erature fall into two classes. One group of workers, using the traditional set 
of microphysics, found spectral differences between 9 e and on the 25% 
level, a range which largely agrees with our findings in view of the different 
microphysics. Other papers claim ratios as large as (ep c ) : (e v ) = 1 : 1.8 or 
even exceeding 1 : 2. We have no explanation for these latter results. At least 
within the framework of our simple power-law models we do not understand 
which parameter could be reasonably adjusted to reach such extreme spectral 
differences. 

In a high-statistics neutrino observation of a future galactic SN one may 
well be able to discover signatures of flavor oscillations. However, when study- 
ing these questions one has to allow for the possibility of very small spectral 
differences, and conversely, for the possibility of large flux differences. This 
situation is almost orthogonal to what has often been assumed in papers 
studying possible oscillation signatures. A realistic assessment of the poten- 
tial of a future galactic SN to disentangle different neutrino mixing scenarios 
should allow for the possibility of very small spectral differences among the 
different flavors of anti-neutrinos. The spectral differences between v e and 
are always much larger, but a large SN neutrino (as opposed to anti- 
neutrino) detector does not exist. 

Previous analyses on oscillation effects in the neutrino flux of a SN not 
only assumed equipartition of the luminosities and a large hierarchy in mean 
energies, but also heavily relied on the absolute magnitude of model pre- 
dictions. Comparing such predictions with results from an observation has 
little predictive power due to uncertainties in the models, e.g., the nuclear 
equation of state. Reliable methods involve only relative flux differences and 
small relative differences in energies. With such rather model independent 
assumptions one can unambiguously identify oscillation effects. 

We proposed a method for identifying earth-matter effects by compar- 
ing the SN neutrino signal of two detectors. One detector records the signal 
from above, the other through the earth. If the energy deposited per unit 
volume in both detectors differs, the earth-matter effect is detected. Sur- 
prisingly, the future IceCube detector at the South Pole, will be able to 
determine the deposited energy with a precision comparable to the discussed 
Hyper-Kamiokande detector. The unique location has the advantage that 
with Super- or Hyper-Kamiokande as a co-detector the likelihood for en- 
countering the required setup (exactly one line of sight through the earth) 
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is 70 %. If the earth- matter effect is detected we can infer the hierarchy of 
neutrino masses to be normal if sin 2 # 13 > 10~ 3 is established by other ex- 
periments. For sin 2 9i 3 < 10~ 3 the earth effect will always be present. In any 
positive signal allows for strong conclusions whereas the reasons for 
a not observing an effect can either be due to smaller flux differences from 
the SN or oscillation parameters. 

Our second proposed mechanism, namely identifying the earth effect at 
a single detector would, together with the first method, enhance the covered 
area of the sky to 85 %. A detector like the proposed Hyper-Kamiokande or a 
future large scintillation detector will record the neutrinos from a galactic SN 
with such a high accuracy that the Fourier transform of the inverse-energy 
spectrum can reveal modulations induced by the earth effect. It would even 
be possible to determine the solar mass splitting at the level of accuracy that 
will be reached by KamLAND. 

These methods have the advantage of being independent of the exact 
original fluxes and spectra. We just assume that a difference between the v e 
and fluxes is present. For the two-detector setup we propose to use the 
time variation of the flux difference as the observable in order to get rid of 
systematic uncertainties in the absolute detector normalization. This time 
dependence is found in all numerical simulations. The physical reason is that 
neutrino emission is first powered by the accretion of infalling material and 
later on by proto-neutron star cooling. 

With this first systematic study of all relevant neutrino-matter inter- 
actions we were able to significantly improve the understanding of flavor- 
dependent SN neutrino emission. Our findings essentially mark a change of 
paradigm for oscillation studies involving SN neutrinos. The mean energies 
of and u c are almost equal whereas the luminosities differ by up to a fac- 
tor of 2, contrary to the previous paradigm. The new picture allows one to 
detect neutrino-oscillation effects, for example with two distant detectors by 
comparing the time variation of the deposited energy, or in a single detector 
by identifying the modulations of the spectrum with the help of a Fourier 
transform. These concepts rely on robust features of SN simulations instead 
of comparing predictions with measurements, like previous methods. The de- 
tection of a SN by itself will be a very valuable astrophysics observation and 
can additionally contribute to particle physics. 
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Abbreviations 



unless specified differently, stands for u T , 9^, and v. 

unless specified differently, stands for and v T 
CM — center of momentum 
LTE — local thermal equilibrium 
Mq — solar mass 
mfp — mean free path 
SN — supernova 
Y c — electron fraction per baryon 
Yl — lepton fraction per baryon 
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Appendix B 
Monte Carlo Code 



B.l General Concept 

Our Monte Carlo code is based on that developed by Janka (1987) where 
a detailed description of the numerical aspects can be found. The code was 
first applied to calculations of neutrino transport in SNe by Janka Sz Hille- 
brandt (1989a,b) and Janka (1991). It uses Monte Carlo methods to follow 
the individual destinies of sample neutrinos (particle "packages" with suit- 
ably attributed weights to represent a number of real neutrinos) on their way 
through the star from the moment of creation or inflow to their absorption or 
escape through the inner or outer boundaries. The considered stellar back- 
ground is assumed to be spherically symmetric and static, and the sample 
neutrinos are characterized by their weight factors and by continuous values 
of energy, radial position, and direction of motion, represented by the cosine 
of the angle relative to the radial direction. The rates of neutrino interactions 
with particles of the stellar medium can be evaluated by taking into account 
Fermion blocking effects according to the local phase-space distributions of 
neutrinos (Janka & Hillebrandt 1989b). 

As background stellar models we use the ones described in Sec. 14.21 They 
are defined by radial profiles of the density p, temperature T, and electron 
fraction Y e , i.e. the number of electrons per baryon. The calculations span the 
range between some inner radius i?i n and outer radius -R ou t- These bound the 
computational domain which is divided into 30 equally spaced radial zones. 
In each zone p, T, and Y e are taken to be constant. R{ n is chosen at such high 
density and temperature that the neutrinos are in LTE in at least the first 
radial zone. R ont is placed in a region where the neutrinos essentially stream 
freely. At i?j n neutrinos are injected isotropically according to LTE. While 
a small net flux across the inner boundary develops, the neutrinos emerging 
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from the star are generated almost exclusively within our computational 
domain. If Ri n is chosen so deep that the neutrinos are in LTE, the assumed 
boundary condition for the flux will not affect the results. 

The stellar medium is assumed to be in thermodynamic equilibrium with 
nuclei being completely disintegrated into free nucleons. Based on p, T, and 
Y c we calculate all the required thermodynamic quantities, notably the num- 
ber densities, chemical potentials, and temperatures of protons, neutrons, 
electrons, positrons, and the relevant neutrinos. Except for runs that include 
weak magnetism, the chemical potentials for and v T are taken to be zero. 
Next we compute the interaction rates in each radial zone for all included pro- 
cesses. In the simulations discussed in the present work, fermion phase-space 
blocking is calculated from the neutrino equilibrium distributions instead of 
the computed phase-space distributions. This simplification saves a lot of 
CPU time because otherwise the rates have to be re-evaluated whenever the 
distribution of neutrinos has changed after a transport time step. The ap- 
proximation is justified because phase-space blocking is most important in 
regions where neutrinos frequently interact and thus are close to LTE. Test 
runs without this approximation show that the results are not affected within 
our numerical accuracy. 

At the start of a Monte Carlo run, 800,000 test neutrinos are randomly 
distributed in the model according to the local equilibrium distributions. 
Each test neutrino represents a certain number of real neutrinos. In this initial 
setup the number of real neutrinos is determined by LTE. Then transport is 
started. The time step is fixed at At = 10~ 7 s; recall that the interaction rates 
do not change. At the beginning of each step neutrino creation takes place. 
The number of test particles that can be created is given by the number 
of neutrinos that were lost through the inner and outer boundaries plus 
those absorbed by the medium. Based on At, the production rates, and the 
fact that the inner boundary radiates neutrinos, we calculate the number of 
neutrinos that are produced in one time step and distribute them among the 
available test neutrinos by attributing suitable weight factors. The sample 
particles are created within the medium or injected at the inner boundary in 
appropriate proportions. 

During a time step the path of each test particle through the stellar 
atmosphere is followed by Monte Carlo sampling. With random numbers and 
stored tables of the total interaction rates we decide whether it flies freely or 
interacts. If it interacts it can scatter or it can be absorbed; in the last case we 
turn to the next particle. For scattering we determine the new momentum 
and position and the time of flight before the interaction took place and 
continue with the process until the time step is used up. Particles leaving 
through the lower or upper boundaries are eliminated from the transport. 
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After a certain number of time steps (typically around 15,000) the neu- 
trino distribution reaches a stationary state and further changes occur only 
due to statistical fluctuations. At that stage we start averaging the output 
quantities over the next 500 time steps. 

B.2 Structure of the Code 

Thermodynamic Properties of the Medium For each radial zone the 
number of baryons is computed by the given medium density. Using the 
given temperature T and Y e together with the number density of baryons 
and the fact that electrons and positrons are in LTE we obtain the chemical 
potential of electrons and positrons by numerically solving Equation (|4.4|) . 
With [l and T the number densities of electrons and positrons are calculated. 
In the same way the properties of protons and neutrons are determined. 

Main Part: The Time-Step Loop The main part of the program con- 
sists of a loop that processes one time step after the other. Each timestep 
corresponds to 10~ 7 s of neutrino transport. Since we use blocking factors of 
neutrinos in LTE as an approximation to the real neutrino occupation, we 
only calculate the interaction rates in the first time step and store them in 
arrays. From these rates the cumulated probabilities are obtained. In the 
time-step loop specific rates are obtained from the arrays by linear interpo- 
lation. 

Injection and Creation of Neutrinos In each time step the number of 
sample neutrinos that were absorbed or left the model through the inner or 
outer boundary during the preceding timestep is now available for injection 
at the inner boundary and emission events in the medium. According to the 
likelihood of each event these sample neutrinos are divided into one group 
that is injected at the inner boundary and various others that correspond to 
the creation mechanisms inside the medium. With the same likelihoods the 
number of real neutrinos that would be produced in the time step is calcu- 
lated and equally distributed among the sample neutrinos as their weights. 
Neutrinos injected at the lower boundary get Ri n as their radial position, a 
random energy according to the LTE distribution in the first cell, and an 
angle according to the distribution of neutrinos freely streaming off a sphere. 
The creation mechanisms inside the medium are the charged current inter- 
actions e~p — > v e ri and e + n — > u e p for v e and 9 e , e + e _ — > vv for all flavors, 
and NN — > NNz/^z^, v e v e — > v^v^ for v^. By Monte Carlo methods we obtain 
the radial position, energy, and angle versus the radial direction. 
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Transport After all new neutrinos are created transport starts. The num- 
ber of neutrinos to be transported equals the number of neutrinos that re- 
mained in the star plus the injected and created neutrinos. The transport 
loop then takes one neutrino after the other through the star until their time 
is used up or they leave the model. Neutrinos that remained inside the model 
from the preceding timestep have the whole time At = lCT 7 s. Since injected 
and created neutrinos appear at some random time their transport time step 
lasts only r At where r is a uniformly distributed random number between 
and 1. Then the time until the first interaction (Tia) is obtained with the help 
of the highest total interaction rate for the energy of the transported neu- 
trino. If Tja lies outside At we store the neutrino properties after it freely 
traveled until the end of At and jump to the next one. If the interaction 
takes place within the available time we have to decide which interaction 
takes place. The interaction rates are also density dependent, i.e. for most 
locations inside the star we underestimated Tja by using the upper limit of 
the total rate. Therefore, when deciding which interaction takes place one of 
the "interactions" leaves the neutrino properties unchanged and in fact is no 
interaction at all. This compensates for the too low estimate Tia locally. 

Before we go to the interaction we first calculate the new position of 
the neutrino at the point of interaction. If it leaves the model through the 
inner or outer boundary on its way to that point it is eliminated and we 
jump to the next one. In case it stays within the area of interest we use 
the interaction rates at the interaction point to decide which process takes 
place. Then according to that process the new energy and angle versus radial 
direction of the outgoing neutrino are determined and stored. If At is not 
used up we determine a new Tja and go on until At is used up. 

After all neutrinos used up their time or left the star we go over to the 
next time step. 
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